Draft version February 2, 2008 

Preprint typeset using I^T^jX style cmulatcapj v. 6/22/04 



INTRINSIC CURVATURE IN THE X-RAY SPECTRA OF BL LACERTAE OBJECTS 
Eric S. Perlman 1,2 , Greg Madejski 3,4 , Markos Georganopoulos 1 ' 5 , Karl Andersson 3,0 , Timothy Daugherty 1 , 



Julian H. Krolik 



Travis Rector 7,8 , John T. Stocke 9 , Anuradha Koratkar 1 
Aller 13 , Hugh Aller 13 , Mark G. Allen 14 

Draft version February 2, 2008 



Stefan Wagner 12 , Margo 



in 
o 
o 

(N 

X> 
<D 
IX, 

m 

(N 
> 
oo 

On 
(N 
(N 
O 

in 
o 

Or 

6 

U 

in 
c3 



ABSTRACT 

We report results from XMM-Newton observations of thirteen X-ray bright BL Lacertae objects, 
selected from the Einstein Slew Survey sample. The survey was designed to look for evidence of 
departures of the X-ray spectra from a simple power law shape (i.e., curvature and/or line features), 
and to find objects worthy of deeper study. Our data are generally well fit by power-law models, with 
three cases having hard (r < 2; dN/dE oc E~ r ) spectra that indicate synchrotron peaks at E £ 5 keV. 
Previous data had suggested a presence of absorption features in the X-ray spectra of some BL Lacs. 
In contrast, none of these spectra show convincing examples of line features, either in absorption or 
emission, suggesting that such features are rare amongst BL Lacs, or, more likely, artifacts caused by 
instrumental effects. We find significant evidence for intrinsic curvature (steepening by dT/d(logE) — 
0.4 ± 0.15) in fourteen of the seventeen X-ray spectra. This cannot be explained satisfactorily via 
excess absorption, since the curvature is essentially constant from 0.5 — 6 keV, an observation which 
is inconsistent with the modest amounts of absorption that would be required. We use the XMM- 
Newton Optical Monitor data with concurrent radio monitoring to derive broadband spectral energy 
distributions and peak frequency estimates. From these we examine models of synchrotron emission 
and model the spectral curvature we see as the result of episodic particle acceleration. 
Subject headings: galaxies: active (galaxies:) BL Lacertae objects: general (galaxies:) BL Lacertae 
objects: individual (1ES0033+595, 1ES0120+340, 1ES0145+138, 1ES0323+022, 
1ES0347-121, 1ES0414+009, 1ES0647+250, 1ES1028+511, 1ES1101-232, 
1ES1133+704 (Mkn 180), 1ES1255+244, 1ES1553+113, 1ES1959+650) X-rays: 
galaxies radiation mechanisms: non-thermal 



1. INTRODUCTION 

Among all active galactic nuclei, BL Lacertae (BL Lac) 
objects are the most dominated by variable, non-thermal 
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emissions. BL Lacs often vary on time scales of days to 
hours at frequencies from optical through the 7-rays (see 
Ulrich, Maraschi & Urry 1997 for a review). Their emis- 
sions are dominated by a broad, featureless continuum, 
believed to originate in a relativistic jet oriented very 
close to our line of sight (see Urry & Padovani 1995 for 
a review). This component is thought to be responsible 
for their variability as well as their high optical polar- 
ization (Jannuzi, Smith & Elston 1994) and radio core- 
dominance (Rector et al. 2000, Rector & Stocke 2001 
and references therein). The spectral energy distribu- 
tions (SEDs) of BL Lacs appear to be dominated by syn- 
chrotron emission at radio to ultraviolet energies (up to 
X-ray energies for X-ray selected objects) and inverse- 
Compton emission at higher energies. 

Previous observations of the X-ray spectra of BL Lacs 
have yielded somewhat confusing results. Worrall & 
Wilkes (1990) compiled the first examination of a rela- 
tively large sample of BL Lacs with Einstein. They found 
power-law shapes, with a wide range of spectral indices. 
This was confirmed by observations with HEAO-1 (Sam- 
bruna et al. 1994), ROSAT (Perlman et al. 1996a, Urry 
et al. 1996, Sambruna et al. 1996), ASCA (Kubo et 
al. 1998) and BeppoSAX (Wolter et al. 1998; Beckmann 
et al. 2002; Padovani et al. 2001, 2004). Those ob- 
servations also found a pattern of harder X-ray spectra 
in objects with spectral peaks in the infrared to opti- 
cal (often called low-energy peaked BL Lacs, or LBLs) 
than in those with spectral peaks in the UV to X-rays 
(often called high-energy peaked BL Lacs, or HBLs). 
The correlation between peak frequency and spectral in- 
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dex is strong, but not monotonic (Padovani, Giommi & 
Fiore 1997, Lamer et al. 1996, Padovani et al. 2001). 
The explanation put forward most often for the complex 
nature of this relationship was that it was due to ob- 
serving the intrinsically curved synchrotron and inverse- 
Compton spectral components, peaked at a variety of 
frequencies, over only a restricted energy range. 

If indeed there is curvature in the X-ray spectra of BL 
Lacs, one might expect to see evidence of this when larger 
bandpasses are examined, particularly in bright objects 
where high signal-to-noise can be attained. Several au- 
thors have found evidence of such curvature (Inoue & 
Takahara 1996; Takahashi et al. 1996; Tavecchio et al. 
1998; Massaro et al. 2004a, b; Giommi et al. 2002) 
among a significant fraction (up to 50%) of HBLs. How- 
ever, those authors did not analyze in detail the alter- 
native possibility of an additional absorbing column, as 
had been assumed by other workers (e.g., Pcrlman et al. 
1996a, Urry et al. 1996, Kubo et al. 1998, Wolter et al. 
1998, Beckmann et al. 2002, Padovani et al. 2001). 

Another issue is whether line features are present in 
the X-ray spectra. Early Einstein grating data for PKS 
2155—304 suggested a soft X-ray deficit (Canizares & 
Kruper 1984), as did later Einstein Solid State Spec- 
trometer data for a few of the brightest objects (Urry, 
Mushotzky & Holt 1986; Madejski et al. 1991). Ex- 
amination of ASCA and BBXRT spectra of some of the 
brightest BL Lacs indicated a similar soft X-ray deficit, 
requiring a recovery of the spectrum towards even lower 
energies, as inferred from ROSAT PSPC data (PKS 
2155-304, Madejski et al. 1992; H1426+428, Sambruna 
et al. 1997; PKS 0548-322, Sambruna & Mushotzky 
1998). These features were all explained by invoking X- 
ray absorption features at 0.5-0.8 keV. However, not all 
bright BL Lacs were found to require such features (e.g., 
Mrk 421, Guainazzi et al. 1999); nor were such features 
found in stacked ROSAT spectra of the fainter EMSS 
BL Lacs (Perlman et al. 1996a). This lack of consensus 
regarding the presence or lack of X-ray spectral curva- 
ture and line features motivated us to use XMM-Newton 
to obtain high signal to noise (S/N) spectra of a signifi- 
cant sample of X-ray bright BL Lac objects with a single 
instrument. 

The paper is organized as follows. In §2, we describe 
the sample, observations and data reduction procedures. 
In §3 we give the results of our X-ray spectral fits, specifi- 
cally concentrating on issues regarding spectral curvature 
and the presence or lack of spectral lines. In §4, we focus 
on models for X-ray spectral curvature and particle ac- 
celeration. Finally, in §5, we conclude with a summary 
and a discussion of the implications of our results. 

2. SAMPLE, DATA AND DATA REDUCTION 

2.1. Sample Design 

Our targets were selected from the Einstein Slew Sur- 
vey sample of BL Lacs (Perlman et al. 1996b). The Slew 
Survey is ideal for this purpose because it is a nearly 
all-sky survey, with mean limiting flux F (0.3-3.5 keV) 
~ 1.4 x 10" 11 erg cm" 2 s" 1 (Elvis et al. 1992). Besides 
being the largest collection of BL Lacs at such high X-ray 
fluxes, the Slew Survey was the first sample to contain 
statistically significant numbers of both HBLs and LBLs. 
We used this property to design a sample of 36 objects, 



which are the X-ray brightest in the LBL and HBL sub- 
classes. Of these, 11 (nearly all HBLs) were on the XMM- 
Newton or Chandra observing lists of other projects and 
so we did not choose to observe them again. We were 
awarded observing time for the 12 X-ray brightest of the 
remaining 24 objects, which unfortunately include only 
the HBL subclass: LBLs were approved only in priority 
C and none were observed, under a separate proposal by 
some members of our team. Table 1 lists the objects we 
observed. We do not discuss here the objects observed 
by other workers, which have already appeared in the 
literature (e.g., Boiler et al. 2001; Watson et al. 2004; 
Blustin, Page & Branduardi 2004; Cagnoni et al. 2004); 
however, where appropriate we make use of their results. 

2.2. Instruments and Observations 

We observed each object with all instruments aboard 
XMM-Newton. The XMM-Newton Observatory (Jansen 
et al. 2001) consists of three coaligned 7.5 m focal 
length X-ray telescopes, focussing X-rays onto the Euro- 
pean Photon Imaging Camera (EPIC) and two Reflection 
Grating Spectrometers (RGS1 and RGS2 respectively). 
The EPIC has three detectors, sensitive in the 0.3-10 keV 
band: two metal-oxide-semiconductor (MOS) and one p- 
n junction (PN) CCD arrays. These instruments pro- 
vide moderate-resolution X-ray imaging (PSF FWHM 
6") and spectroscopy (E/5E « 20 — 50, depending on 
energy) . The two RGS instruments provide spectroscopy 
at higher resolution, with E/SE w 300. 

Because the main goal was to address the sample prop- 
erties of the BL Lac class, rather than to obtain very 
deep observations of individual objects, our integration 
times were short, ~ 5 ks on each instrument. Table 1 
lists the EPIC PN and MOS good on-source times and 
count rates. Figure 1 shows two examples of our EPIC 
observations, fit with a single-power-law model (§3). As 
can be seen, the S/N of the EPIC data were quite high 
(20-50 below 1 keV; lower at high energies). The RGS 
spectra (§3.3) had much more modest S/N, due to the 
lower effective area of the gratings. In all but one ob- 
servation, the EPIC was used in small-window mode in 
order to minimize the effects of pileup. On most objects, 
the THIN filter was used in both the PN and MOS; how- 
ever, for the brighter objects we used the MEDIUM filter 
to avoid contamination from optical light due to the ob- 
ject (none of the targets has other X-ray bright objects 
within the field). 

Eight XMM-Newton observations were significantly af- 
fected by proton flares. Of these, six were impacted badly 
enough that one or more of the X-ray instruments had 
to be shut off because of the high background, and in 
two, all instruments were shut off during the observa- 
tion. For the six observations where X-ray instruments 
were shut off completely due to high background, the ob- 
servation was automatically attempted again (see Table 1 
for details). As a result, four objects were observed twice 
and one object was observed three times. We inspected 
all data from XMM-Newton, regardless of whether the 
observation was affected by flares. We generally found 
that even for the observations impacted by proton flar- 
ing, some of the X-ray data were usable as long as the 
instruments were not shut off. 

XMM-Newton also has a coaligned 30 cm optical/UV 
telescope, the Optical Monitor (OM). The OM is sensi- 
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tive between A = 1000 — 6000 A and has a typical an- 
gular resolution of ~ 2". OM data were gathered for 
all objects, but in only one band due to observing re- 
strictions in AOl that did not allow multi-band obser- 
vations within < 2 hours. All but one OM dataset was 
gathered using the UVW2 filter, which is sensitive be- 
tween 2000 - 3000 A. Data were taken in exposures of 
800 s in order to allow for the possibility of finding short- 
timescale optical variability, but no such variations were 
found. We include in Table 2 the near-UV magnitude of 
each source, as derived from the OM data. 

We also monitored the radio fluxes of all but one of 
these objects using the University of Michigan Radio Ob- 
servatory (UMRAO) 26-m telescope (Aller et al. 1985). 
We observed each object at one or more of three frequen- 
cies: 4.5, 8.0 and 14.5 GHz, but some were below its sen- 
sitivity limitations (typically 50 mJy at 8.0 GHz). These 
data are also given in Table 2; two example datasets are 
shown in Figure 2. 

2.3. Data Reduction Procedures 

All source and background extraction, as well as exam- 
ination of the lightcurves, was done in the XMM-Newton 
Science Analysis System (SAS) v5.4.1. We performed 
a standard reduction of the events list for the PN, MOS 
and RGS data, which involves the subtraction of hot and 
dead pixels, removal of events due to electronic noise, and 
the correction of event energies for charge transfer losses. 
The PN files were also filtered to include only single, dou- 
ble, and single+double events (PATTERN < 4) while the 
MOS included all single to quadruple events (PATTERN 
< 12). Photons were extracted in all instruments using a 
45" aperture (which, for the PN and MOS, was circular 
in shape). In cases where pile- up was significant, an an- 
nular extraction region was used to filter out the pixels 
affected most severely. 

We extracted lightcurves for every dataset, and filtered 
out all high-background time intervals. To do so, we de- 
termined the count rate in the 10 - 15 keV band for the 
entire PN and MOS detectors. Since the source count 
rate in this band is quite small, intervals with more than 
0.3 ct s~ x for the PN and 1.5 ct s _1 for MOS were re- 
jected. Once this was done, the datasets were searched 
for evidence of source variability on < 5 ks timescales 
(using the nominal 0.3 - 10 keV band); no such variabil- 
ity was found during any of these observations. Spectra 
were extracted, as discussed in the next section. These 
spectra were limited to the energy range 0.2-10 keV for 
the MOS and 0.15-15 keV for the PN; however, for the 
detailed spectral fits, we used a more restricted energy 
range, as given below, to account for possible residual 
calibration uncertainties (see Andersson & Madejski 2004 
for detailed discsusion). 

For the MOS data, the background region used was 
the same size as the source region in a place where there 
are no visible sources (typically off axis). For the PN 
data, the background region was once again the same 
size as the source region at a point where the RAWY 
values were the same as the source. Because of the high 
fluxes of these sources and the use of the small-window 
mode (as well as the relatively stringent criteria to re- 
ject data segments affected by flares), we did not extract 
background maps (see e.g., Read & Ponman 2003; note 
that the background maps vary little within the portion 



of the detectors used for small- window mode). It should, 
however, be noted that because of the use of the small- 
window mode the background regions are of necessity 
close to the source itself, and thus it is possible that some 
small amount of source flux could be subtracted. This 
effect causes mainly a small error in the absolute flux 
of the sources, and thus does not affect our analysis of 
possible curvature or line features. We did, however, at- 
tempt to minimize this effect by placing the background 
region as far away from the source as possible within the 
small window. 

For the OM data, we performed simple aperture pho- 
tometry, with a 6" circular aperture, using an annulus 
for background subtraction. The size of the annulus was 
typically 20"; however, in the case of a few objects, a 
nearby companion necessitated some adjustment. We 
then used the zero-points published in the XMM-Newton 
Observer's Handbook (Issue 2.1) to calculate absolute op- 
tical fluxes. The UMRAO data were reduced according 
to the procedure outlined in Aller ct al. (1985). The 
fluxes and magnitudes from UMRAO and the OM at the 
time of observation are summarized in Table 2. 

3. X-RAY SPECTRAL FITS 

As the sources were well centered in the XMM observ- 
ing region, we used standard response (RMF) matrices 15 , 
picking the one closest in time to each observation. The 
Ancillary Response Files (ARFs) were generated using 
the arfgen function of SAS 5.4.1. To assure the valid- 
ity of Gaussian statistics, we grouped the data, combin- 
ing the instrumental channels such that each new bin 
would have at least 40 counts. The exceptions are 1ES 
0145+138 (OBSID 0094383401) and 1ES 1255+244 (OB- 
SID 0094383201), where because of the small number of 
counts we rebinned to 10 counts/bin for the MOS data 
and 15 counts/bin for the PN data. All the PN data 
were fit in the 1.1-10.0 keV range, while all the MOS 
data were fit in the 0.3-10.0 keV range, except for those 
of the faintest objects which we capped at 7.0 keV due 
to low count rates at the highest energies. 

All X-ray spectral modeling was done in XSPEC vll.0. 
Continuum models were fit only to the PN and MOS 
data and then applied to the RGS data. We used both 
the EPIC and RGS data to search for absorption and/or 
emission lines (§3.3). The PN and MOS data were fit si- 
multaneously as well as individually in order to spot any 
calibration problems and/or instrument specific prob- 
lems that might have occurred during the observation. 
Where multiple observations of an object were obtained 
(Table 1, §2.2), each observation was reduced and ana- 
lyzed separately, as we expected that the flux and X-ray 
spectrum would change between those observations, as 
was indeed the case. 

3.1. Spectral Modeling in the broad X-ray band 

Initially, two models were fit to each EPIC spectrum: 

1. A single power law, 

dN/dE = ke- a W N "<°°>e- a W NH ^ l+ ^ E~ T '. 

2. A logarithmic parabola, 

dN/dE = fce-^^Yff.Gaie-^E^H.^tU+z) E {-v+pLo a {E)) _ 

15 These files are available from 

ftp: / / xmm.vilspa.esa.es / pub/ccf/constitucnts/cxtras /responses / 
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In each case, we attempted fits with two different treat- 
ments of the absorbing column: with Nu,int = and 
with Nn,int (assumed to be at the redshift of the source) 
allowed to vary freely. We took values of Nn,Gal horn 
Stark et al. (1992), although it is probable that the ac- 
tual Galactic absorbing column for each object might be 
slightly different due to the low resolution of the Stark 
et al. survey. In all cases, the opacity associated with 
the intervening column density was modeled by a stan- 
dard Morrison & McCammonn (1983) absorption model. 
The motivation to try both the power law model, as well 
as the less-used logarithmic parabola model (described 
in Giommi et al. 2002 and references therein), was that 
both simple power-law and continuously curving spec- 
tral shapes can reasonably be expected depending on 
the model adopted for synchrotron aging and acceler- 
ation (Leahy 1991, Massaro 2002). Our spectral fits for 
the single power-law model are summarized in Table 3. 
In Figure 1, we show as examples spectra of two objects 
fit with the single power-law model. 

As can be seen, the spectral indices we found range 
from r = 1.7 — 2.9 (corresponding to energy spectral in- 
dices a = 0.7 — 1.9), with twelve of seventeen being in the 
range T = 2.1 — 2.7. These spectral indices are similar 
to the findings of past observations of HBLs by ROSAT 
(Pcrlman al. 1996a, Padovani et al. 1997), ASCA (Kubo 
et al. 1998) and BeppoSAX (Wolter et al. 1998, Beck- 
mann et al. 2001). Three spectra were found to be flat 
(r < 2). These objects are likely to have synchrotron 
peaks at energies higher than ~ 5 keV (see also §4). 

In the three objects where repeated observations were 
made, we see significant variability both in flux and spec- 
tral shape (Table 3). In all of those objects, the high- 
est flux state observed also has the hardest spectrum, 
in agreement with the known spectral variability prop- 
erties of BL Lac objects (e.g., Ulrich et al. 1997 and 
references therein). However, the correlation between a 
harder spectrum and higher X-ray flux is not one-to-one, 
as indicated by the three observations of 1ES1959+650. 

3.2. X-ray Spectral Curvature 

In fourteen of seventeen observations, a better fit (sig- 
nificant at the > 99% level according to F-test results; 
see Table 3) was obtained by allowing for the possibility 
of spectral curvature (over and above the default, power- 
law plus Galactic Nh model). This curvature can either 
be intrinsic, or the result of additional absorption due to 
material either within the BL Lac's host galaxy or our 
own. As indicated above, we first investigated the possi- 
bility that additional absorption was present; the result 
of these procedures is given in Table 3. As can be seen, 
in 13 of 14 observations the column required is Si 20% 
of the Galactic figure. We also investigated the alternate 
possibility of intrinsic curvature, using the logarithmic 
parabola model as well as an alternate model, detailed 
later, where we simply fitted power-laws independently 
in four sub-bands, and obtained results of similar signifi- 
cance. The results of this procedure are given in Table 4 
and discussed below. Importantly, it is impossible to dis- 
criminate a priori between the two models due to their 
similar mathematical forms (see above). It is, however, 
possible to use other information and perform other tests. 

We first inspected the implied columns to check the 
reasonability of the idea that the spectral curvature was 



the result of additional absorbing material. We believe 
that additional absorbing material is less likely for three 
reasons. First, the typical BL Lac's host galaxy is a 
bright elliptical, where columns this high would not nor- 
mally be expected except in extreme cases (Goudfrooij 
et al. 1994). Also, in 1ES1959+650, one of the three 
objects where repeated observations were done, the im- 
plied Nh appears to vary between epochs. This requires 
the absorbing material to have been in a region smaller 
than ~ 1 light-month in size. This would imply densities 
n £ 10 3 — 10 4 cm~ 3 , which is outwardly not unreason- 
able for the inner regions of the AGN. But the implied 
density would have to climb with decreasing timescale, 
as t -3 / 2 , meaning that unreasonable column densities 
would be reached if Nh appeared to vary on timescales 
only 1-2 orders of magnitude smaller, which are not un- 
common for flux variability among BL Lacs (Ulrich et al. 
1997). Finally, there is no indication from optical or UV 
spectroscopy that BL Lac objects would have any signif- 
icant neutral material in the line of sight to the nucleus 
(e.g., Perlman et al. 1996; Rector et al. 2000; Rector 
& Stocke 2001; Kinney et al. 1991; Lanzetta, Turnshck 
& Sandoval 1993; Penton & Shull 1996). This material, 
in principle, can be partially ionized, providing no opac- 
ity in the optical and UV, but absorbing via individual 
edges in the soft X-ray band. This can be ruled out as a 
general property on the basis of an absence of any indi- 
vidual spectral features in the RGS data (§3.3), although 
we cannot exclude it in the spectra of fainter objects. 

This is persuasive but not conclusive evidence that 
these spectra have intrinsic X-ray spectral curvature. To 
test this idea further, we split the XMM-Newton data 
into four sub-bands: 0.5-1 keV, 1-2 keV, 2-4 keV and 4- 
10 keV, and fit each sub-band individually with a power- 
law model. For this exercise we assumed an Nh value 
fixed at Galactic, and used the MOS data only below 
2 keV and both the PN and MOS data at higher en- 
ergies. This was done because of the known issues in 
cross-calibrating PN data with those from the MOS at 
lower energies, especially near 1 keV (see also the first 
paragraph of §3, above). The results of this procedure 
are given in Table 4, and four examples of these fits are 
shown in Figure 3. We note that the overall \ 2 values 
for these fits are comparable to those achieved by varying 
Nh; however it is impossible to use an F-test to compare 
them because of the differing ways in which data was se- 
lected for this procedure. As can be seen, the typical 
X-ray curvature seen in our sample is nearly constant, 
even up to the highest energy band. We have char- 
acterized this curvature by looking at the dT / dilogE) 
for each object. This information is included in Ta- 
ble 4. As can be seen, our data are consistent with a 
(dT / d(\ogE)) w 0.4±0.15 in fourteen of seventeen obser- 
vations (the same ones for which spectral curvature was 
indicated as mentioned above). This is not consistent 
with the observed columns in these objects, which range 
from 3 — 20 x 10 20 cm~ 2 . At such moderate columns, it 
is typical for the spectrum above about 1 keV to be rel- 
atively insensitive to the absorbing column. This is yet 
another argument (albeit indirect) against the presence 
of additional absorbing material explaining the observed 
spectral curvature in these objects. 

As already noted, we also attempted to model intrin- 
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sic curvature by fitting a logarithmic parabola form. The 
information for these fits is also given in Table 4. The 
logarithmic parabola is characterized by both a spectral 
index and a curvature parameter, /?. The spectral form 
is somewhat different mathematically from either of the 
other forms; however, it is essentially similar to the ap- 
plication of a constant dT/d(logE), with the added ad- 
vantage of allowing a consistent fit in the entire band. In 
theory, the value of the curvature parameter (5 should be 
equivalent to the value of (dT /dilogE)) found using sev- 
eral sub-bands. Inspection of Table 4 shows that these 
are indeed similar, albeit with small variations in each 
case. The mean curvature required for the objects where 
this is indicated by the F-tests is essentially identical, in 
fact: (/?) = 0.37 versus (d/T/d(logE)) = 0.38. 

Giommi et al. (2002) also found good evidence for 
intrinsic X-ray spectral curvature (which they modeled 
by the same logarithmic parabola form) in about half 
of the HBLs observed by BeppoSAX. The value of (ft) 
they found is consistent with our findings. However, the 
proportion of curved spectra in Giommi et al. (2002) is 
considerably lower than we see. The difference cannot 
be explained by our small sample, as twelve of thirteen 
objects observed by BeppoSAX were also observed by us 
with XMM-Newton, and of those, for six the best repre- 
sentation of their BeppoSAX spectrum was a logarithmic 
parabola plus fixed absorption (although n.b., Giommi et 
al. 2002 did not fit a variable absorption model for any 
object). We believe it is much more likely that the reason 
for the larger fraction found here lies in the higher sig- 
nal to noise and resolution of these spectra compared to 
those from BeppoSAX. This is supported by the fact that 
if one examines the variable Njj fits, the 6 objects where 
Giommi et al. (2002) required a logarithmic parabola 
are best fitted by significantly higher absorbing columns 
than those which Giommi et al. (2002) was able to fit 
adequately with a single power law. 

Our conclusion from the above is that we believe the 
most appropriate model for the spectral curvature we 
see is intrinsic, rather than due to additional absorption. 
Given the similar values of x 2 , our primary motivation 
for concluding this is the indirect lines of evidence pre- 
sented above regarding the implausibility of additional 
absorbing columns under these circumstances. We will 
proceed on this basis in the upcoming discussion. 

To close this discussion of intrinsically curved spectra, 
it is useful to compare the fits and shape of the spec- 
trum achieved for a given object under a simple power 
law model (with Galactic Nh) with those implied by the 
logarithmic power law form. We show this comparison in 
Figure 4, for one of the spectra of 1ES1959+650, specif- 
ically OBSID 0094383301. There comparison is shown 
in two forms: at top, we show the comparison in terms 
of the detected spectrum (folded into the instrumental 
response), while at bottom it is shown in terms of emit- 
ted energy. The plots at top enable us to see the real 
differences in the quality of the fit: the simple power 
law model (Figure 4, top left) clearly overpredicts the 
observed spectrum at E > 1 keV, and slightly undcr- 
predicts the lower-noise MOS data at 1-1.5 keV. But at 
higher energies it is the PN data which have the greater 
ability to distinguish between the two models, and the 
simple power law model overpredicts the data at energies 
higher than about 5 keV. As can be seen, the logarith- 



mic power law model (Figure 4, top right) fits the data 
much better. The plots at bottom enable us to see the 
form of the curvature: note that above ~ 1 keV essen- 
tially no curvature would be seen, if the spectrum fol- 
lowed a simple power law model (Figure 4, bottom left). 
However, under the logarithmic parabola model (Figure 
4, bottom right), we see that there is in fact significant 
curvature seen at higher energies, albeit of a mild, grad- 
ual form. This is presumably why this curvature was 
not seen by earlier satellites, which suffered from either 
a much smaller spectral range (for all satellites except 
BeppoSAX) and/or much lower signal to noise (all other 
satellites). 

3.3. Spectral Line Features in BL Lacs? 

We inspected each spectrum thoroughly for spectral 
features. We did not find convincing evidence for line 
features in any of these spectra. To quantify this, both 
the EPIC and RGS data for the 12 highest signal-to- 
noise spectra were searched for absorption features in 
XSPEC (we did not search datasets where either instru- 
ments were turned off and/or the statistics were inad- 
equate). The best- fit broadband model (Table 4) was 
used as the seed model in each case, and the RGS data 
were grouped in bins with a minimum of 20 photons 
apiece (following our practice for the EPIC data; §3.1). 
We then artificially set a line feature of 5 eV physical 
width, and covering fraction of 0.5, and then stepped it 
through the range 0.5 to 4.5 keV in increments of 0.01 
keV using STEPPAR, and assessed the significance of 
features using the A\ 2 (A 3cr feature has A% 2 = —9.21, 
while a 4cr feature has A\ 2 = —18.42). Higher ener- 
gies were not searched due to the small effective area 
and poor statistics. Once features were found, we then 
fixed the energy and allowed XSPEC to converge on a 
best-fit width. This procedure found six narrow fea- 
tures, as shown in Table 5, in four of the datasets, but 
in most (1ES 0120+340 OBSID 0094382101, 1ES0323 
OBSID 0094382501, 1ES 0414+009 OBSID 0094383101, 
1ES0647+250 OBSID 0094380901, 1ES1133+704 OB- 
SID0094170101, 1ES1255+244 OBSID 0094383001, 1ES 
1553+113 OBSID 0094380801, 1ES1959+650 OBSID 
0094383501), no features were found. 

Are these features real? As can be seen, all but one are 
between 3-4 cr, with one barely over 4er. For Gaussian 
statistics, P(3cr) = 0.0027, while P(4cr) = 0.000318. As 
can be seen from Table 3, the typical EPIC dataset in our 
sample is ~ 300 channels, as is the typical RGS dataset 
(n.b., the number of channels is not significantly greater 
due to the relatively poor statistics resulting from the 
short exposure times). Thus the result of 5 features at 
3cr and 1 feature at > 4cr in ~ 3600 spectral channels 
is consistent with what one would expect from statistical 
noise (respectively 10 + 3 features at > 3er and 1 + 1 
feature at > 4er). Moreover, as noted in Table 5, none 
of the features found in any of the spectra are present 
in the spectra from all the EPIC and RGS instruments. 
Thus, we believe all are the result of random noise and 
therefore not real spectral lines. The limits shown in 
Table 5 can be taken as representative of the detection 
limits of our data, which range from ~ 5 up to ~ 50 eV 
depending on count rate. 

This conclusion differs from findings based on data 
from earlier instruments. We are confident of our result 
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based on the relatively high signal-to-noise of the XMM- 
Newton spectra. We note, however, that narrow, small 
equivalent width absorption lines such as those inferred 
(at intermediate redshifts) from Chandra grating obser- 
vations of BL Lac objects by Nicastro et al. (2002; see 
also Cagnoni et al. 2004), or those due to ionized oxygen 
in the hot halo of our Galaxy (e.g., Nicastro et al. 2002, 
McKernan, Yaqoob & Reynolds 2004), are not excluded 
by the data in-hand. The determination of presence or 
absence of such features require much higher S/N data 
than are available in our sample. 

It is worth noting that a similar absence of strong 
spectral features was also found in another, independent 
study of the XMM-Newton X-ray spectra of four of the 
five BL Lacs where BBXRT and ASCA spectra appeared 
to show these line features, specifically H1219+301, Mkn 
501, H1426+428 and H0548-322 (Blustin et al. 2004). 
Those authors reported that if any of the previously 
found absorption features were indeed real, they had 
to represent a transient phenomenon. They found this 
particular conclusion unlikely, ruling it out at 93% con- 
fidence based on the assumption that the (multiple) 
XMM-Newton observations of those objects represented 
random time instances. The objects in the Blustin et al. 
(2004) sample are not discussed here, and moreover, our 
integration times are smaller so we have used primarily 
the EPIC data and not the RGS. However, our results 
are in the same vein. With the addition of our data, it 
appears unlikely that the absorption features found with 
ASCA and BBXRT were transient phenomena. 

The most likely explanation of the previously reported 
line features is a combination of mis-calibration of the 
previous instruments and the use of overly simplistic 
spectral models. As an example, a deficit of low energy 
counts in the ASCA data has been previously noted in 
reference to the data for NGC 5548 (Iwasawa, Fabian, 
and Nandra 1999); the simultaneous ASCA and ROSAT 
PSPC observation implied significantly higher flux be- 
low 0.5 keV than any reasonable spectral extrapolation 
of the ASCA SIS data would allow. Another reason is 
that a spectral curvature of exactly the type discussed in 
this paper would, at low resolution and lower signal-to- 
noise, mimic an absorption feature when the underlying 
spectrum was assumed to be a simple power law. How- 
ever, the confirmation that an absorption feature was 
measured would require data below the energy of such a 
putative feature, and such data were not always available 
(as was the case for the Einstein SSS as well as SSS + 
MPC data). We therefore suggest that the claims made 
by earlier workers were a byproduct of the particular 
simple power law model used, while the effect was due 
to gradual spectral curvature and/or instrumental effects 
rather than absorption features from material along the 
line of sight. 

4. MODELING THE X-RAY AND BROADBAND 
SPECTRAL CHARACTERISTICS 

We have constructed broadband spectral energy dis- 
tributions (SEDs) for each object using the EPIC data, 
in conjunction with the data from the OM and the UM- 
RAO flux nearest in time to the XMM observation. For 
objects which were too faint in either the optical or radio 
at the time of observation to register a positive detection 
on either the OM or UMRAO, we have used POSS plates 



to obtain an average optical flux and radio fluxes from 
the literature (Perlman et al. 1996b and sources therein) 
for the SED. We fit to these data a simple parabolic 
model, in order to obtain a peak frequency, v pe ak- We 
did not assume a standard synchrotron spectrum for this 
procedure because the small amount of data we have do 
not adequately constrain these models (see e.g., Leahy 
1991, Pacholczyk 1970). For a few objects, where the 
initial fit was not good, we had to modify a single non- 
simultaneous point to find a convincing value of v pea k ■ 

In Table 6 we summarize the broadband spectral data 
and Up ea k for every object in our sample. As can be 
seen, our objects range between log v pea k ~ 14.5 — 18, as 
expected for HBL-type and intermediate objects. A few 
objects, noted in Table 6, have very flat X-ray spectral 
index and/or significant evidence of v pea k > 10 18 Hz. 
Our procedure yielded v pea k values considerably in excess 
of 10 18 Hz for all these objects (see also §3.1), but in 
Table 6 we list them as having v pea k = 10 18 Hz because 
our data are not able to test for a higher peak frequency 
given their frequency space coverage. 

We then proceeded to test for correlations be- 
tween broadband spectral properties and both V and 
dT/d(logE). We then tested each for correlations using 
standard statistical tests (Spearman's rank correlation, 
p, and Pearson's correlation coefficient, r). These plots 
and their physical implications are the main subject of 
this section. 

In Figure 5 we show plots of L versus a ro (left), 
a ox (middle) and v pea k- As can be seen, our data 
show significant anti-correlations between T and v pea k 
(r = -0.741, p = -0.781, P = 0.18%) and between P 
and a ox (r = 0.609, p = 0.634, P = 1.1%) but no signifi- 
cant correlation between T and and a ro (r = 0.036, p — 
— 0.05,P = 84.0%). This constellation of results is con- 
sistent with the a x - u pea k relation of Padovani et al. 
(1997), when one considers that the range of v pea k val- 
ues we cover is only ~ 3 decades. 

4.1. Statistical Particle Acceleration? 

The observed spectra exhibit curvature over a fre- 
quency range from 0.5 keV to 10 keV, in the sense 
of a gradual steepening with increasing energy. As- 
suming that the observed emission is synchrotron ra- 
diation, this implies that the electron distribution re- 
sponsible for this emission is not the typical power law 
predicted by particle acceleration schemes (for a recent 
review see Gallant 2002). Instead, the electron spec- 
trum has to gradually curve downwards in a way that 
produces the observed spectra. Given that X-ray syn- 
chrotron emission requires in situ particle acceleration 
due to the short radiative lifetimes of the emitting par- 
ticles (T sync = 1.2 x 10 3 Bc 3/2 £keV 2(5 ~ 1/2 s ; R-ybicki & 
Lightman 1979), the ultimate physics behind curved X- 
ray spectra must be connected on a very basic level with 
the nature of particle acceleration in these systems. Here 
and in the next subsection we discuss possible solutions 
to this problem. 

In Massaro et al. (2004b) a log-parabolic spectral cur- 
vature was fit to the X-ray spectrum and broadband SED 
of Mkn 421 in several flux states. Those authors mod- 
eled the log-parabolic curvature via a statistical particle 
acceleration process under which a particle i with energy 
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7i in a given region has a finite probability (< 1) of being 
accelerated, where the probability is given by pi = g/rf, 
where g, q are positive constants. They show in their 
paper that this leads analytically to a log-parabolic cur- 
vature form (see their §6 and equations (10)-(18)), and 
predicts a linear correlation between the X-ray spec- 
tral index, T x , and dT/d(logE). In Figure 6, we show 
plots of Vpeak versus dT /d(logE) (top) as well as T x ver- 
sus dT /d(logE) (bottom). As can be seen, neither of 
these plots shows a significant correlation (y pea ^ versus 
dT/d{logE) : r = 0.615, p = 0.393, P = 11.6% ; 1^ ver- 
sus dT/d{logE) : r = -0.24,/) = 0.014, P = 95.3%). 
The results are similar for tests of J^peak 

and T versus (3. 

Thus our data do not support the model of Massaro et 
al. (2004b). One possible reason why this is so is that 
Massaro et al. (2004b) model particle acceleration only, 
and do not include losses. Moreover, the statistical par- 
ticle acceleration is still assumed to be time-invariant at 
all locations in the jet and hence the issue of variability 
is not addressed. 

4.2. Spectral Curvature: a Signature of Episodic 
Particle Acceleration? 

The results of §4.1 prompt us to look for another expla- 
nation that takes into account the known observational 
properties of jets and particularly blazars. Another mo- 
tivation for this is the recent result of Pcrlman & Wilson 
(2005) that standard, continuous injection models of par- 
ticle acceleration consistently overpredict the X-ray flux 
of the M87 jet by large factors (up to ~ 100 at 1 keV, 
varying with energy roughly as E 0A ). The interpreta- 
tion advanced in that paper was one of a position- and 
energy-variable (but not time-variable) filling factor for 
particle acceleration. The Perlman & Wilson work, how- 
ever, does not explore the issue of spectral curvature, 
as the statistics in most components of the M87 jet are 
inadequate for this purpose; nor does it explore issues re- 
lated to variability, which are known to be of paramount 
importance for BL Lac objects. 

In this subsection we address these issues directly, by 
demonstrating that curved electron particle spectra can 
be produced if the particle acceleration is episodic. Im- 
portantly, this model assumes only a time-variable par- 
ticle acceleration which in an integrated sense should be 
indistinguishable from a a spatial average. Unlike the 
case of M87, here we have no spatial information regard- 
ing the distribution of jet X-ray emission and/or particle 
acceleration. Instead, we know that variability is an im- 
portant characteristic of blazar jets and thus we optimize 
our model for this case. We explore these subjects more 
deeply in a later paper (Georganopoulos & Perlman, in 
preparation) that fully outlines a model considered in 
brief form below. 

Consider a zone, possibly a shock, where electrons with 
Lorentz factor 7 are injected at a rate of Q electrons 
per second and accelerated to higher energies. Following 
Kirk, Rieger, & Mastichiadis (1998), the kinetic equation 
describing particle acceleration is 
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= Q*(7-7o), (1) 



where n(7,i) is the electron energy distribution, l/t acc 
is the particle acceleration rate, 1 /t esc is the escape rate 



of particles from the acceleration region, and 

0=1 — U (2) 
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with ot the Thomson cross-section and U the total en- 
ergy density of the magnetic field and ambient photons. 
Assuming that the injection started at time t = 0, the 
electron energy distribution after time t will be a power 
law with electron index s = 1 + t acc jt esc up to a Lorentz 
factor 
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where "/ max — l//3t acc is the maximum energy electrons 
reach asymptotically if the particle acceleration mecha- 
nism operates for t > \n[(j max - ja)/jo}t a cc- The time 
t required for electrons to be accelerated up to Lorentz 
factor 71 is 



, , , lljlmax ~ 7o) 

t(ji)-t acc ln— — — . 
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If we assume that particle acceleration is characterized 
by a typical timescale T, then the electron distribution 
at any time < t < T is a power law that cuts off at a 
Lorentz factor 7(f), reaching its maximum Lorentz factor 
71 (T) at t = T. Particle acceleration ceases to operate 
at t = T, and the total time r over which electrons of 
Lorentz factor 7 are provided by the particle acceleration 
mechanism is 



T — t = t. 



, 7l(7max -7) 
acc nl . . . 

7(7maz - 7l) 



(5) 



The time r is also the time during which an electron of 
Lorentz factor 7 radiates, as long as the radiative cooling 
timescale t rac i = l//?7 = t acc j max /j satisfies the condi- 
tion t rad <C t. In other words, r is the total time the 
accelerated electron distribution is non-zero at electron 
Lorentz factor 7; note that, as is apparent from Eq. (5), 
the condition t rac i -C t can always be satisfied for 71 
sufficiently close to j max - If the light crossing time of 
the emission region or the integration time of our obser- 
vations is greater than T, then we effectively observe a 
time averaged electron distribution {n(j)) which is the 
product of a power law multiplied by a logarithmic term 
representing the time interval this power law is available 
at each energy: 



(71.(7)} oc 7 s In 



71 (7n 
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7 < 7i- 



(6) 



lilmax -7l)' 

Using the the 5- function approximation (i.e., an electron 
of Lorentz factor 7 produces synchrotron photons only 
at the critical synchrotron energy) to calculate the syn- 
chrotron emission, we obtain an observed photon spec- 
trum 



dN 



ocS-( s+1 )/ 2 ln^ 
dE E^{EUl 



1/2, 



E<E U (7) 



where Ei and E max are the synchrotron energies pro- 
duced by electrons of Lorentz factor 71 and "f ma ,x re- 
spectively. In Figure 7 we show the photon spectrum 
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and photon index as a function of energy, for the fol- 
lowing parameter choices: s = 3.0, E max = 100 keV and 
Ei = 50kcV. As can be seen, the spectrum steepens grad- 
ually, with T increasing at a fairly constant rate up to ~ 
few keV, and then dramatically at E > 10 KeV. If our 
interpretation of the steepening is correct, we expect to 
see an increasing curvature in future observations above 
the 10 keV upper limit of XMM-Newton, as predicted 
in equation (7) (see Figure 7). One might also expect 
that in other objects, e.g., LBL, this cutoff might occur 
at lower energies. In LBL, however, we do not observe 
such steep values of T, but rather, flat X-ray spectra are 
seen (F < 2; see Padovani et al. 2001, 2004), which 
are typically interpreted as being dominated by inverse- 
Comtpon emission (Padovani et al. 1997, Lamer et al. 
1996). In practice, it may be difficult to observe this 
steepening because as the curvature of the synchrotron 
component increases, the spectrum is essentially cutting 
off and at some point the observed spectrum will begin 
to be dominated by the onset of the Compton compo- 
nent. What one might expect to see in practice, then, 
would be concave spectra (i.e., increasing T) up until the 
two components cross, with an inflection point at their 
crossing point and then a constant value of T (< 2) at 
higher energies. One final thing that should be noted 
here is that above 10 keV, almost all the HBL spectra 
that currently exist are for objects that are currently in 
a flaring state, where we would not expect to see signif- 
icant steepening because the particle acceleration is still 
at or near its maximum. Also, there are no currently 
available instruments that would be sufficiently sensitive 
at E > 10 keV to search for further steepening in the 
spectrum of HBLs. 

The two energies Ex and E max appearing in Equation 
(7) can be derived through spectral fitting. However, 
their values alone are offering little insight to the phys- 
ical conditions in the source, because they depend on 
a series of parameters (B, t acc , 70, T). As we discuss 
in §5, additional information from frequency-dependent 
variability can be used to constrain the physical state of 
the source. 

5. CONCLUSIONS 

We have presented 17 high-signal-to-noise XMM- 
Newton X-ray spectra for 13 BL Lac objects. These 
spectra do not show any spectral lines, contrary to ex- 
pectations from some previous observations. We believe 
the reason for this discrepancy is due to a combination 
of overly simplistic spectral models and/or insufficiently 
precise calibration of previous instruments. As a result 
the intrinsic X-ray spectra of BL Lac objects appear to 
be quite featureless, giving no sign as to the physics of 
any region within the central engine, other than the jet 
itself. This is quite different from all other non-blazar 
type sources, where some line features, in either emission 
or absorption, are seen. This can be interpreted in the 
context of unified scheme models (e.g., Urry & Padovani 
1995), under which the jet emission is Doppler boosted 
because it is seen at small angles to our line of sight. It 
is, however, difficult to address the issue of the necessary 
viewing angles to produce the degree of featurelessness 
seen in these spectra, particularly given that at present 
our knowledge concerning the X-ray spectra of the par- 
ent population of FR 1 radio galaxies is limited. Fu- 



ture work will provide much information on this subject, 
but currently the best information available (based on 
luminosity function estimates) suggests viewing angles 
9 ~ 10 — 30 and Lorentz factors r ~ 5. This would yield 
Doppler boosts to the jet's apparent luminosity of hun- 
dreds to thousands, as detailed in e.g., Urry & Padovani 
(1995). 

Our data are best fit by curved spectra which grow 
steeper in an approximately logarithmic fashion across 
the XMM-Newton band. From an observational point 
of view, a logarithmic curvature is not a new concept, 
having been first introduced by Landau et al. (1986) 
for radio-optical observations of blazars, and then re- 
vived for X-ray spectra by Giommi et al. (2002). Our 
analysis is, however, the first to demonstrate that in- 
trinsic curvature is to be preferred over the alternate 
model of additional absorbing material along the line of 
sight. From the theoretical point of view, logarithmically 
curved spectra require explanation, as most of the stan- 
dard models of synchrotron emission produce spectra of 
a constant power-law slope (see, e.g., Leahy 1991 for a 
review). Massaro et al. (2004b) suggested statistical ac- 
celeration, whereby a given particle in the distribution 
has a probability of being accelerated, as one possible 
model for this type of spectral shape. The predictions of 
that model are not borne out by our analysis, however, 
and a statistical model also does not take into account 
the variable nature of BL Lacs. It was with this moti- 
vation that we considered an episodic model for particle 
acceleration, which explicitly considers a modification to 
the kinetic equation and thus includes both losses and a 
model of variability. We showed that this model fits our 
data reasonably well and predicts continued steepening 
at higher energies, until the onset of the IC component. 

Episodic particle acceleration may be related to the 
broadband variability observed in TeV sources like 
1959+650 (Krawczynski et al. 2004), where both the 
power and peak frequency of the synchrotron spectrum in 
X-rays exhibit an increase by more than a factor of ~ 10, 
while the optical flux remains practically constant. This 
can be easily accommodated in our scheme through an 
increase of the characteristic time T that particle acceler- 
ation operates. In this case 71 and E\ will both increase, 
resulting in a more powerful X-ray synchrotron spectrum 
that peaks at higher energies. The situation will be quite 
different at lower frequencies; the curvature of the spec- 
trum produced by the logarithmic term will be present 
down to photon energies for which the radiative cooling 
times t r ad becomes larger than the characteristic parti- 
cle acceleration time T. Below this energy the electron 
energy distribution and the corresponding synchrotron 
spectrum will be a simple power law. At these lower en- 
ergies an increase in T can only affect the amplitude of 
the observed spectrum and even that only to the extent 
that such an increase results in an increased fraction of 
time that particle acceleration is operating (often called 
the duty cycle). If, for example, T describes both the 
time acceleration is on and off, then as T increases the 
average number of particles injected remains constant, 
resulting in a constant flux at lower frequencies, while 
the spectrum changes dramatically at high frequencies 
due to the increase of both 71 and Ex . It is very encour- 
aging that our observations of 1ES1959+650, one to three 
months after flare observed by Krawczynski et al. (2004), 
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show significant spectral curvature. This indicates that 
the characteristic timescales on which episodic acceler- 
ation might operate is ^ 1 month, which is consistent 
with the findings of long-look campaigns (Perlman et al. 
1999, Tanihata et al. 2001). Such a timescale is also con- 
sistent with Doppler-boosted versions of the flare found 
recently in the jet of M87 (Harris et al. 2003, Perlman 
et al. 2003), requiring only values S ~ 10 (as opposed to 
the more modest 5=1 — 2 required for M87). Indeed, 
the fact that we see logarithmic curvature in 14 of 17 of 
these spectra indicates that the typical ratio of r/T ;$ 1; 
i.e., the amount of time in which flares are seen is not 
appreciably smaller than the amount of time in which 
more gentle or no variability is seen. This is not inconsis- 
tent with the data from the ASCA long-look campaigns 
(Tanihata et al. 2001). 
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Fig. 1— Two examples of the XMM EPIC spectra. At top, Mkn 180 (OBSID 0094170101) and at bottom, 1ES1959+650 (OBSID 
0094383301). In both panels, the EPN data are the top line and the data from the MOS detectors are below. The model shown represents 
the power law plus variable Njj fit. See Sections 2 and 3 for discussion. 
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Fig. 2. — Radio light curves for 1ES 1553+113 (top) and 1ES 1028+511 (bottom), obtained at the University of Michigan Radio 
Astronomy Observatory. Data at 14.5, 8.0 and 4.8 GHz are denoted by crosses, circles and triangles respectively. Daily averages of the 
data are shown, with the bars denoting la error estimates. As seen, both sources vary considerably during the radio observations. Sec 
Section 2.1-2.2 for discussion. 
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Fig. 3. — Four examples of fitting the spectra of BL Lacs in our sample to smaller spectral ranges. All four of these objects show steeper 
spectra at higher energies, with dT/d(logE) Ri 0.4. See Section 3.3 for discussion. 
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Fig. 4.— Comparison of model fits for the XMM spectrum of 1ES1959+650 (OBSID 0094383301). In the panel at left, the data were 
fitted with a power-law plus Galactic N(H) model, while at right, the data were fitted with the log-parabola plus Galactic N(H) model. 
Note that the fit in the left panel is adequate between 1.5-5 keV, but cannot account for the curvature seen both at energies < 1.5 keV and 
> 5 keV. At top, the comparison is shown in terms of the detected spectrum, with each plot showing received spectra plus folded model 
for the MOS 1 & 2 (bottom lines) and PN (top lines), as well as residual. At bottom, we show the unfolded energy spectrum. These plots 
show that curvature is present in this spectrum up to the highest energies seen by XMM. 
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Fig. 5. — Plots of the X-ray spectral slope T, versus broadband spectral index a ro (top left) and a ro (top right), as well as the SED peak 
frequency u pea ^ (bottom left). For easy cross-referencing with the Tables, we note in these plots the identity of each point. See Section 4.1 
for discussion. 
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Fig. 6.— Plots of dT/d{logE) versus T (top) and /^p ea fc (bottom). For easy cross-referencing with the Tables, we note in these plots the 
identity of each point. See Section 4.1 for discussion. 
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Fig. 7. — Upper panel: The photon spectrum (solid line) resulting from an episodic injection of electrons in the particle acceleration 
process. Here s = 3.0, E max = 100 KeV, and E\ = 50 KeV. Note that E\ is the maximum emitted energy the source is reaching in this 
episodic mode of injection, while E max is the emitted energy the source would reach under steady injection conditions. For comparison, 
a power law emission corresponding to an electron distribution with s = 3.0 is also shown (broken line). Lower panel: the photon index 
T that corresponds to the emitted photon spectrum. Note the gradual steepening of the spectrum, manifested by an increase of T with 
energy. The dotted lines represent the energy range 0.5 — 10 keV covered by XMM. See Section 4.2 for discussion. 



TABLE 1 
Log of XMM Observations 



1ES Name 


OBSID 


Date 


MOS gt 


Filter 


Frame 


Ct Rate 


PNgt 


Filter 


Frame 


Ct Rate 








(sec) 






(sec) 


(sec) 






(Ct/s) 


0033+595 


0094381301 1 


2003-02-01 


1500 


Thin/Med 


PartialW2 


1.17 


4100 


Thinl 


Small 


3.43 


0120+340 


0094382101 


2002-01-05 


5600 


Medium 


PartialW2 


2.43 


3500 


Medium 


Small 


8.48 


0145+138 


0094383401 


2002-07-18 


4200 


Thinl 


PartialW2 


0.07 


3200 


Thinl 


Small 


0.25 




AAA A ooocm 

UUy4ooz5Ul 


zUUz-Uz-U5 


5200 


Medium 


Partial W 2 


0.81 


3200 


Medium 


Small 


2.81 


no a t 1 01 
Uo4 f — 1Z1 


nnn /iqoi i m 
L)Uy4oollUl 


onno no oc 


3000 


Thinl 


Partial W 2 


3.32 


2200 


1 mill 


omall 


11.6 


ri/i 1 a i nnn 
U414-|-UUy 


nno/i qqq i m 
UUy4ooolUl 


onno ns ok 

ZUUZ-Uo-ZO 


o t UU 


i nmi 


rartiai/ run 


QIC 
O. lo 


ODUU 


Tallin 1 

i nmi 


omaii 


1Z. 1 


nfi/iv i oe;n 
U04 f -j-ZOK) 


uuy4oouyui ' 


onno riQ 
ZUUZ-Uo-Zo 




Medium 


Partial Wz 


o i n 
z. 1U 


i nnn 


Medium 


omall 


*7 a n 






9009 03 9^ 


N /A 


N /A 


N /A 


N /A 


N /A 


M 1 1 1 m 
IvlCtll Hill 


Small 


N /A 




uuy'iooioui 


9nni 1 ^ 
zuui-uo-io 


N /A 


N /A 


N /A 


N /A 




Medium 


Olildll 


o.oy 




0094382701 


2001-11-26 


5600 


Medium 


PartialW2 


4.30 


3500 


Medium 


Small 


16.2 


1101-232 


0094380601 1 


2001-05-29 


2200 


Medium 


PartialW2 


4.70 


4300 


Medium 


Small 


17.6 


1133+704 


0094170301 1 


2001-04-12 


N/A 


Thinl 


Full 


N/A 


N/A 


Thinl 


Full 


N/A 




0094170101 2 


2001-04-12 


3300 


Thinl 


Full 


4.45 (2.02) 3 


6900 


Thinl 


Full 


12.5 (5.81) 3 


1255+244 


0094383001 1 


2002-12-12 


5800 


Medium 


PartialW2 


1.30 


4000 


Medium 


Small 


4.40 




0094383201 


2002-06-26 


662 


Medium 


PartialW2 


1.35 


435 


Medium 


Small 


4.88 


1553+113 


0094380801 


2001-09-06 


5000 


Medium 


PartialW2 


13.0 (4.31) 3 


3500 


Medium 


Small 


50.5 


1959+650 


0094380201 1 


2002-11-23 


300 


Medium 


PartialW2 


21.2 (5.29) 3 


N/A 


Medium 


Small 


N/A 




0094383301 1 


2003-01-16 


2200 


Medium 


PartialW2 


10.5 (2.08) 3 


1200 


Medium 


Small 


39.5 




0094373501 


2003-02-09 


3700 


Medium 


PartialW2 


12.5 (3.80) 3 


2000 


Medium 


Small 


43.3 



x This observation was affected by proton flaring 

2 Observation was attempted again in the same orbit 

3 Count Rates are raw count rates in the effective exposure interval, while values in parentheses are rates after pile-up correction. 
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TABLE 2 
Supporting observations 



1ES Name OBSID OM Band F(UV) -F max (rad) 









(10 10 erg cm z s 1 A i ) 


(Jy) 


0033+595 


0094381301 


UVW2x2 


< 3.3 1 


0.3 


0120+340 


0094382101 


UVW2x5 


8.9 ± 0.9 


0.1 


0145+138 


0094383401 


U VW2 x 4 


< 3.3 1 


0.2 


0323+022 


0094382501 


UVW2x5 


14.4 ± i.o 


0.3 


f)347— 1 21 


0094381 1 01 


VJ V VYiAl 


7.6 ± 2.0 


2 


041 zi4-nno 


ooqa^&^i ni 


T TVW9 v \ 

U V VV Zr A O 


io.o in u. ± 


1 

u. 1 


nfi/i i i oc;n 
U04 / -f-ZOU 


nno/i Qcnoni 
UUy4ooUyUl 


U V W Z X 1 


< 0.0 


n 1 
U. 1 




0094382901 


UVW2 x 5 


< 3.3 1 


0.1 


1028+511 


0094381801 


U VW2 x 5 


18.2 ± 1.6 


0.2 




0094382701 


UVW2 x 4 


21.9 ± 0.7 


0.2 


1101-232 


0094380601 


U VW2 x 4 


10 ± 3 


0.3 


1133+704 


0094170301 


N/A 


N/A 


2 




0094170101 


Ux4 


272 ±8 


N/A 3 


1255+244 


0094383001 


U VW2 x 5 


3.2 + 1.5 


0.2 




0094383201 


U VW2 x 4 


< 3.3 1 


0.2 


1553+113 


0094380801 


U VW2 x 5 


185 + 5 


0.8 


1959+650 


0094380201 


U VW2 x 5 


35.0 + 2.0 


0.6 




0094383301 


U VW2 x 5 


25.3 ± 2.0 


0.6 




0094373501 


U VW2 x 5 


28.4 + 2.0 


0.6 



1 Object not detected in the OM observation; value is a 2cr upper limit 

2 Below detection threshold in all UMRAO observations 

3 No UMRAO observations made as this object was originally part of another program 



TABLE 3 
BL Lac Spectral Fit Parameters 



1ES Name 


Obsld 




Fit parameters, N(H) Fixed 






Fit parameters, N(H) Free 








N H ,Gal 


T F2-10keV 


X 2 /dof. (dof.) 


N H ,tot 


r F 2 _ 10kcV X 2 /dof. (dof.) 


F-tcst 






10 20 cm" 2 


erg cm~ 2 s _1 




10 20 cm" 2 


erg cm 2 s 1 


Probability 



00^+^9^ 


0094^81 ^01 F 


42.4 


9 9+0.02 
z ' z -0.03 


1.2 




10" 


11 


1 4.9(9%%} 


f?n q+3.8 


9 47+0 0^ 


1.1 x 10" 


n 


1 01 (9R7\ 


<? in -10 


01 20+340 


0094382101 


5.14 


2 16+ ' 01 


9.7 


x 


10" 


12 


1 70('527 > ) 


9.9±0.6 


2 34+0- 02 
z -° -0.03 


8.9x10" 


12 


1 32('526 > ) 


< io- 10 


01 4^-Ul ^8 

U14iJT lOO 


00QzL383zL01 


5.10 


9 07 + 19 


2.5 




10" 


13 




99 1+8.0 


9 no+0.38 


1.6x10" 


13 


1.24(87) 


3.2 x 10 - 4 


0393-1-099 


ooqa3&9^oi 


R 7 A 
o. ( 4 


9 c-4+0.02 


9 




"1 o - 


12 


1 0^91 3^ 


11 9±1 2 


9 fi7+0 OR 


1 9 x 10" 


12 


n 0^91 o\ 


1 3 x 10" ^ 


0347—121 


0094381 101 


3.64 


i 77 +0.02 
1 - ' ' -0.01 


1.8 


x 


10" 


11 


1 ISfSSO 1 ) 


5.0±0.6 


1 82+0 03 


1.8x10" 


11 


1 lfi('52Q s ) 


2.6 x 10 -3 


0414+009 


0094383101 


10.3 


9 47 +0.02 

z - 4 '-o.oi 


9.3 


X 


io- 


12 


1.20(499) 


13.1±0.6 


2.57±0.03 


8.9x10- 


12 


1.09(498) 


< IO" 10 


0o4l +250 


0094380901 


12.8 


9 ^9 _l_ n n9 
z.bz ± 0.02 


5.0 


X 


10 


12 


1.30(2 to ; 


19 -°±i'.i 


2.8o±0.05 


5.0x10" 


12 


1.05(2 //) 


^ in — 10 
< 1U 


1028+511 


0094381801 F 


1.16 


-0.03 


1.2 


X 


io- 


11 


1.20(107) 


3.77±1.1 


2.61±0.08 


1.0x10" 


11 


1.04(106) 


< IO" 10 




0094382701 


1.16 


2.28 ±0.01 


1.2 


X 


io- 


11 


1.41(600) 


5.2±0.4 


2.46±0.02 


9.7x10- 


12 


0.98(599) 


< IO" 10 


1101-232 


0094380601 F 


5.76 


9 90+O.O2 

z - zo -o.oi 


2.3 


X 


io- 


11 


1.37(682) 


10.9±0.6 


2 40+ ' 02 
z - 4u -0.03 


2.2x10- 


11 


1.05(681) 


< IO" 10 


1133+704 


0094170101 


1.42 


2.40 ± 0.02 


9.8 


X 


io- 


12 


1.21(361) 


1.8±0.7 


2.42±0.04 


9.7x10- 


12 


1.20(360) 


8.4 xlO- 2 


1255+244 


0094383001 F 


1.26 


1 99+ 01 


5.1 


X 


io- 


12 


1.18(380) 


2.9± 2 2 9 6 


1.83±0.12 


4.8x10- 


12 


1.02(379) 


< IO" 10 




0094383201 


1.26 


1 76+ 07 
x - ,D -0.06 


7.5 


X 


io- 


12 


0.97(134) 


4.6±0.8 


9 -i i +0.04 
z - 11 -0.03 


4.8x10- 


12 


0.97(133) 


N/A 


1553+513 


0094380801 


3.67 


2.38 ±0.01 


3.5 


X 


io- 


11 


1.26(815) 


5.7±0.4 


2 46+ ' 01 
z - 4D -0.02 


3.4x10- 


11 


1.17(814) 


< io- 10 


1959+650 


0094380201 F ' a 


10.1 


1.70 ±0.05 


2.9 


X 


io- 


10 


1.44(58) 


13 3+ 31 


l-80-o.ii 


2.8x10- 


10 


1.41(57) 


0.275 




0094383301 F 


10.1 


2 10+ - 02 


7.7 


X 


io- 


11 


1.51(514) 


17.9±1.0 


9 9Q+0.02 
^•^ a -0.03 


7.2x10- 


11 


1.08(513) 


< io- 10 




0094383501 


10.1 


2.02 ± 0.02 


6.9 


X 


io- 


11 


1.64(863) 


16.9±0.6 


9 1O+0.02 

Als -o.oi 


6.5x10- 


11 


1.13(862) 


< IO" 10 



Note. — MOS data used in 0.5 - 10.0 keV range PN in 1.1 - 10.0 keV range. 

a No PN data available, fit parameters are for MOS data only. 
F Most of this observation is contaminated by proton flaring. 



TABLE 4 

Parameters for BL Lac Spectral Curvature 



4-band T-Fitting 12 Logarithmic Parabola™ 



Lhio iNamc 


UDsld 


J- 0.5-1.0fceV 


1 1.0-2.0fceV 


1 2.0-4.0fceV 


1 4.0-10. OkeV 


XZ/aoi. (aoi.j 


al / (lylogrj ) 


1 


a 
P 


^,0 JA^f ( A^f \ 

XZ/dol. (aoi.j 


0033+595 


0094381301 


i-49ti;° 9 I 


2 03+ ' 18 


2 S4+ ' 12 
Z "J^_0.12 


q q o-hO . 22 

Z " 3J — 0.22 


1.04(187) 


0.58 


1 71 01 

lm ' -0.09 


-o.67i8:l 2 


1.10 (287) 


0120+340 


0094382101 


i.94i8;i 


2 14 +0.08 
z - ±4 -0.08 


2 19+° 1 
z - ±M -0.09 


q on +0.22 
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"'Nh assumed fixed at Galactic (see Table 3) for all fits. 

b X 2 /dof. figures were derived by separately adding the x 2 and dof figures and only then dividing. 
c Four-band not performed because only PN data of low S/N were available 



TABLE 5 

Discrete features found in RGS Spectra 



1ES Name" 


OBSID 


Energy (keV) 


Ax 2 


Width (cV) 


Comments 


1ES0033+595 
1ES1101-232 
1ES1101-232 
1ES1255+244 
1ES1255+244 
1ES1959+650 


0094381301 
0094380601 
0094380601 
0094383201 
0094383201 
0094383301 


1.07 
0.67 
1.36 
1.29 
2.91 
1.92 


-12.02 
-18.56 
-16.62 
-11.98 
-8.91 
-11.77 


45 +20 

40 -24 
7 +4 
'-3 
22 ±9 
30 ± 16 

8±« 
< 18 


MOS 2 + 1 RGS only 
RGS 2 only 
RGS 2 only 

MOS 1, MOS 2 but not in RGS 

Single MOS 1&2 channel, instrumental? 

Two RGS 1 channels, not RGS 2 or MOS 


a This table lists only spectra 


where discrete features 


were found. The spectra 


where the search was 


performed but no features were found 


are listed in §3.3 
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MULTIWAVEBAND PARAMETERS 
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0094381801 
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1101-232 


0094380601 


0.45 


0.67 
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1133+704 


0094170301 


0.26 


1.38 
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0.48 


0.72 


18.00 




0094373501 


0.48 


0.79 


18.00 



Comments 



Converges poorly; radio flux had to be lowered artificially 



Confusing radio source? Catalog flux 25% of UMRAO value 
a x < 1 - use log v peak =18. 

Catalog optical flux 15x brighter than OM data 
Catalog optical flux 15x brighter than OM data 



Catalog opt flux 25x brighter than OM; increased opt flux 2x 
a x < 1 - use log i/ peak = 18. 

a x < 1 - use log v peak = 18. 

Used log y pea k = 18 even though object not in flare. 
Used log Vp e ak = 18 even though object not in flare. 
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ABSTRACT 

We report results from XMM-Newton observations of thirteen X-ray bright BL Lacertae objects, 
selected from the Einstein Slew Survey sample. The survey was designed to look for evidence of 
departures of the X-ray spectra from a simple power law shape (i.e., curvature and/or line features), 
and to find objects worthy of deeper study. Our data are generally well fit by power-law models, with 
three cases having hard (T < 2; dN/dE oc E~ T ) spectra that indicate synchrotron peaks at E 5 keV. 
Previous data had suggested a presence of absorption features in the X-ray spectra of some BL Lacs. 
In contrast, none of these spectra show convincing examples of line features, either in absorption or 
emission, suggesting that such features are rare amongst BL Lacs, or, more likely, artifacts caused by 
instrumental effects. We find significant evidence for intrinsic curvature (steepening by dT/d(\ogE) = 
0.4 ± 0.15) in fourteen of the seventeen X-ray spectra. This cannot be explained satisfactorily via 
excess absorption, since the curvature is essentially constant from 0.5 — 6 keV, an observation which 
is inconsistent with the modest amounts of absorption that would be required. We use the XMM- 
Newton Optical Monitor data with concurrent radio monitoring to derive broadband spectral energy 
distributions and peak frequency estimates. From these we examine models of synchrotron emission 
and model the spectral curvature we see as the result of episodic particle acceleration. 
Subject headings: galaxies: active (galaxies:) BL Lacertae objects: general (galaxies:) BL Lacertae 
objects: individual (1ES0033+595, 1ES0120+340, 1ES0145+138, 1ES0323+022, 
1ES0347-121, 1ES0414+009, 1ES0647+250, 1ES1028+511, 1ES1101-232, 
1ES1133+704 (Mkn 180), 1ES1255+244, 1ES1553+113, 1ES1959+650) X-rays: 
galaxies radiation mechanisms: non-thermal 
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Among all active galactic nuclei, BL Lacertae (BL Lac) 
objects are the most dominated by variable, non-thermal 
emissions. BL Lacs often vary on time scales of days to 
hours at frequencies from optical through the 7-rays (see 
Ulrich, Maraschi & Urry 1997 for a review). Their emis- 
sions are dominated by a broad, featureless continuum, 
believed to originate in a relativistic jet oriented very 
close to our line of sight (see Urry & Padovani 1995 for 
a review). This component is thought to be responsible 
for their variability as well as their high optical polar- 
ization (Jannuzi, Smith & Elston 1994) and radio core- 
dominance (Rector et al. 2000, Rector & Stocke 2001 
and references therein). The spectral energy distribu- 
tions (SEDs) of BL Lacs appear to be dominated by syn- 
chrotron emission at radio to ultraviolet energies (up to 
X-ray energies for X-ray selected objects) and inverse- 
Compton emission at higher energies. 

Previous observations of the X-ray spectra of BL Lacs 
have yielded somewhat confusing results. Worrall & 
Wilkes (1990) compiled the first examination of a rela- 
tively large sample of BL Lacs with Einstein. They found 
power-law shapes, with a wide range of spectral indices. 
This was confirmed by observations with HEAO-1 (Sam- 
bruna et al. 1994), ROSAT (Perlman et al. 1996a, Urry 
et al. 1996, Sambruna et al. 1996), ASCA (Kubo et 
al. 1998) and BeppoSAX (Wolter et al. 1998; Beckmann 
et al. 2002; Padovani et al. 2001, 2004). Those ob- 
servations also found a pattern of harder X-ray spectra 
in objects with spectral peaks in the infrared to opti- 
cal (often called low-energy peaked BL Lacs, or LBLs) 
than in those with spectral peaks in the UV to X-rays 



2 



(often called high-energy peaked BL Lacs, or HBLs). 
The correlation between peak frequency and spectral in- 
dex is strong, but not monotonic (Padovani, Giommi & 
Fiore 1997, Lamer et al. 1996, Padovani et al. 2001). 
The explanation put forward most often for the complex 
nature of this relationship was that it was due to ob- 
serving the intrinsically curved synchrotron and inverse- 
Compton spectral components, peaked at a variety of 
frequencies, over only a restricted energy range. 

If indeed there is curvature in the X-ray spectra of BL 
Lacs, one might expect to see evidence of this when larger 
bandpasses are examined, particularly in bright objects 
where high signal-to-noise can be attained. Several au- 
thors have found evidence of such curvature (Inoue & 
Takahara 1996; Takahashi et al. 1996; Tavecchio et al. 
1998; Massaro et al. 2004a, b; Giommi et al. 2002) 
among a significant fraction (up to 50%) of HBLs. How- 
ever, those authors did not analyze in detail the alter- 
native possibility of an additional absorbing column, as 
had been assumed by other workers (e.g., Perlman et al. 
1996a, Urry et al. 1996, Kubo et al. 1998, Wolter et al. 
1998, Bcckmann et al. 2002, Padovani et al. 2001). 

Another issue is whether line features are present in 
the X-ray spectra. Early Einstein grating data for PKS 
2155—304 suggested a soft X-ray deficit (Canizares & 
Kruper 1984), as did later Einstein Solid State Spec- 
trometer data for a few of the brightest objects (Urry, 
Mushotzky & Holt 1986; Madejski et al. 1991). Ex- 
amination of ASCA and BBXRT spectra of some of the 
brightest BL Lacs indicated a similar soft X-ray deficit, 
requiring a recovery of the spectrum towards even lower 
energies, as inferred from ROSAT PSPC data (PKS 
2155-304, Madejski et al. 1992; H1426+428, Sambruna 
et al. 1997; PKS 0548-322, Sambruna & Mushotzky 
1998). These features were all explained by invoking X- 
ray absorption features at 0.5-0.8 keV. However, not all 
bright BL Lacs were found to require such features (e.g., 
Mrk 421, Guainazzi et al. 1999); nor were such features 
found in stacked ROSAT spectra of the fainter EMSS 
BL Lacs (Perlman et al. 1996a). This lack of consensus 
regarding the presence or lack of X-ray spectral curva- 
ture and line features motivated us to use XMM-Newton 
to obtain high signal to noise (S/N) spectra of a signifi- 
cant sample of X-ray bright BL Lac objects with a single 
instrument. 

The paper is organized as follows. In §2, we describe 
the sample, observations and data reduction procedures. 
In §3 we give the results of our X-ray spectral fits, specifi- 
cally concentrating on issues regarding spectral curvature 
and the presence or lack of spectral lines. In §4, we focus 
on models for X-ray spectral curvature and particle ac- 
celeration. Finally, in §5, we conclude with a summary 
and a discussion of the implications of our results. 

2. SAMPLE, DATA AND DATA REDUCTION 

2.1. Sample Design 

Our targets were selected from the Einstein Slew Sur- 
vey sample of BL Lacs (Perlman et al. 1996b). The Slew 
Survey is ideal for this purpose because it is a nearly 
all-sky survey, with mean limiting flux F (0.3-3.5 keV) 
~ 1.4 x 10" 11 erg cm" 2 s" 1 (Elvis et al. 1992). Besides 
being the largest collection of BL Lacs at such high X-ray 
fluxes, the Slew Survey was the first sample to contain 



statistically significant numbers of both HBLs and LBLs. 
We used this property to design a sample of 36 objects, 
which are the X-ray brightest in the LBL and HBL sub- 
classes. Of these, 11 (nearly all HBLs) were on the XMM- 
Newton or Chandra observing lists of other projects and 
so we did not choose to observe them again. We were 
awarded observing time for the 12 X-ray brightest of the 
remaining 24 objects, which unfortunately include only 
the HBL subclass: LBLs were approved only in priority 
C and none were observed, under a separate proposal by 
some members of our team. Table 1 lists the objects we 
observed. We do not discuss here the objects observed 
by other workers, which have already appeared in the 
literature (e.g., Boiler et al. 2001; Watson et al. 2004; 
Blustin, Page & Branduardi 2004; Cagnoni et al. 2004); 
however, where appropriate we make use of their results. 

2.2. Instruments and Observations 

We observed each object with all instruments aboard 
XMM-Newton. The XMM-Newton Observatory (Jansen 
et al. 2001) consists of three coaligned 7.5 m focal 
length X-ray telescopes, focussing X-rays onto the Euro- 
pean Photon Imaging Camera (EPIC) and two Reflection 
Grating Spectrometers (RGS1 and RGS2 respectively). 
The EPIC has three detectors, sensitive in the 0.3-10 keV 
band: two metal-oxide-semiconductor (MOS) and one p- 
n junction (PN) CCD arrays. These instruments pro- 
vide moderate-resolution X-ray imaging (PSF FWHM 
6") and spectroscopy (E/SE « 20 — 50, depending on 
energy) . The two RGS instruments provide spectroscopy 
at higher resolution, with E/5E « 300. 

Because the main goal was to address the sample prop- 
erties of the BL Lac class, rather than to obtain very 
deep observations of individual objects, our integration 
times were short, ~ 5 ks on each instrument. Table 1 
lists the EPIC PN and MOS good on-source times and 
count rates. Figure 1 shows two examples of our EPIC 
observations, fit with a single-power-law model (§3). As 
can be seen, the S/N of the EPIC data were quite high 
(20-50 below 1 keV; lower at high energies). The RGS 
spectra (§3.3) had much more modest S/N, due to the 
lower effective area of the gratings. In all but one ob- 
servation, the EPIC was used in small-window mode in 
order to minimize the effects of pileup. On most objects, 
the THIN filter was used in both the PN and MOS; how- 
ever, for the brighter objects we used the MEDIUM filter 
to avoid contamination from optical light due to the ob- 
ject (none of the targets has other X-ray bright objects 
within the field). 

Eight XMM-Newton observations were significantly af- 
fected by proton flares. Of these, six were impacted badly 
enough that one or more of the X-ray instruments had 
to be shut off because of the high background, and in 
two, all instruments were shut off during the observa- 
tion. For the six observations where X-ray instruments 
were shut off completely due to high background, the ob- 
servation was automatically attempted again (see Table 1 
for details). As a result, four objects were observed twice 
and one object was observed three times. We inspected 
all data from XMM-Newton, regardless of whether the 
observation was affected by flares. We generally found 
that even for the observations impacted by proton flar- 
ing, some of the X-ray data were usable as long as the 
instruments were not shut off. 
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XMM-Newton also has a coaligned 30 cm optical/UV 
telescope, the Optical Monitor (OM). The OM is sensi- 
tive between A = 1000 — 6000 A and has a typical an- 
gular resolution of ~ 2". OM data were gathered for 
all objects, but in only one band due to observing re- 
strictions in AOl that did not allow multi-band obser- 
vations within < 2 hours. All but one OM dataset was 
gathered using the UVW2 filter, which is sensitive be- 
tween 2000 - 3000 A. Data were taken in exposures of 
800 s in order to allow for the possibility of finding short- 
timescale optical variability, but no such variations were 
found. We include in Table 2 the near-UV magnitude of 
each source, as derived from the OM data. 

We also monitored the radio fluxes of all but one of 
these objects using the University of Michigan Radio Ob- 
servatory (UMRAO) 26-m telescope (Aller et al. 1985). 
We observed each object at one or more of three frequen- 
cies: 4.5, 8.0 and 14.5 GHz, but some were below its sen- 
sitivity limitations (typically 50 mJy at 8.0 GHz). These 
data are also given in Table 2; two example datasets are 
shown in Figure 2. 

2.3. Data Reduction Procedures 

All source and background extraction, as well as exam- 
ination of the lightcurves, was done in the XMM-Newton 
Science Analysis System (SAS) v5.4.1. We performed 
a standard reduction of the events list for the PN, MOS 
and RGS data, which involves the subtraction of hot and 
dead pixels, removal of events due to electronic noise, and 
the correction of event energies for charge transfer losses. 
The PN files were also filtered to include only single, dou- 
ble, and single+double events (PATTERN < 4) while the 
MOS included all single to quadruple events (PATTERN 
< 12). Photons were extracted in all instruments using a 
45" aperture (which, for the PN and MOS, was circular 
in shape). In cases where pile- up was significant, an an- 
nular extraction region was used to filter out the pixels 
affected most severely. 

We extracted lightcurves for every dataset, and filtered 
out all high-background time intervals. To do so, we de- 
termined the count rate in the 10 - 15 keV band for the 
entire PN and MOS detectors. Since the source count 
rate in this band is quite small, intervals with more than 
0.3 ct s^ 1 for the PN and 1.5 ct s^ 1 for MOS were re- 
jected. Once this was done, the datasets were searched 
for evidence of source variability on < 5 ks timescales 
(using the nominal 0.3 - 10 keV band); no such variabil- 
ity was found during any of these observations. Spectra 
were extracted, as discussed in the next section. These 
spectra were limited to the energy range 0.2-10 keV for 
the MOS and 0.15-15 keV for the PN; however, for the 
detailed spectral fits, we used a more restricted energy 
range, as given below, to account for possible residual 
calibration uncertainties (see Andersson & Madejski 2004 
for detailed discsusion). 

For the MOS data, the background region used was 
the same size as the source region in a place where there 
are no visible sources (typically off axis). For the PN 
data, the background region was once again the same 
size as the source region at a point where the RAWY 
values were the same as the source. Because of the high 
fluxes of these sources and the use of the small-window 
mode (as well as the relatively stringent criteria to re- 
ject data segments affected by flares), we did not extract 



background maps (see e.g., Read & Ponman 2003; note 
that the background maps vary little within the portion 
of the detectors used for small- window mode). It should, 
however, be noted that because of the use of the small- 
window mode the background regions are of necessity 
close to the source itself, and thus it is possible that some 
small amount of source flux could be subtracted. This 
effect causes mainly a small error in the absolute flux 
of the sources, and thus does not affect our analysis of 
possible curvature or line features. We did, however, at- 
tempt to minimize this effect by placing the background 
region as far away from the source as possible within the 
small window. 

For the OM data, we performed simple aperture pho- 
tometry, with a 6" circular aperture, using an annulus 
for background subtraction. The size of the annulus was 
typically 20"; however, in the case of a few objects, a 
nearby companion necessitated some adjustment. We 
then used the zero-points published in the XMM-Newton 
Observer's Handbook (Issue 2.1) to calculate absolute op- 
tical fluxes. The UMRAO data were reduced according 
to the procedure outlined in Aller et al. (1985). The 
fluxes and magnitudes from UMRAO and the OM at the 
time of observation are summarized in Table 2. 

3. X-RAY SPECTRAL FITS 

As the sources were well centered in the XMM observ- 
ing region, we used standard response (RMF) matrices 16 , 
picking the one closest in time to each observation. The 
Ancillary Response Files (ARFs) were generated using 
the arfgen function of SAS 5.4.1. To assure the valid- 
ity of Gaussian statistics, we grouped the data, combin- 
ing the instrumental channels such that each new bin 
would have at least 40 counts. The exceptions are 1ES 
0145+138 (OBSID 0094383401) and 1ES 1255+244 (OB- 
SID 0094383201), where because of the small number of 
counts we rebinned to 10 counts/bin for the MOS data 
and 15 counts/bin for the PN data. All the PN data 
were fit in the 1.1-10.0 keV range, while all the MOS 
data were fit in the 0.3-10.0 keV range, except for those 
of the faintest objects which we capped at 7.0 keV due 
to low count rates at the highest energies. 

All X-ray spectral modeling was done in XSPEC vll.0. 
Continuum models were fit only to the PN and MOS 
data and then applied to the RGS data. We used both 
the EPIC and RGS data to search for absorption and/or 
emission lines (§3.3). The PN and MOS data were fit si- 
multaneously as well as individually in order to spot any 
calibration problems and/or instrument specific prob- 
lems that might have occurred during the observation. 
Where multiple observations of an object were obtained 
(Table 1, §2.2), each observation was reduced and ana- 
lyzed separately, as we expected that the flux and X-ray 
spectrum would change between those observations, as 
was indeed the case. 

3.1. Spectral Modeling in the broad X-ray band 

Initially, two models were fit to each EPIC spectrum: 

1. A single power law, 

dN/dE = ke-< 7 ^ NH ^'e-^ E *> N »^ 1+z *> E~ r '. 

2. A logarithmic parabola, 

16 These files are available from 

ftp: / / xmm.vilspa.esa.es / pub/ccf/constituents/extras /responses / 
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dN/dE = fc e -< T ( B ) Ar ff,Gai e -<T(B)A r ff,i„t(l+^) E (-r+f3Log(E)) _ 

In each case, we attempted fits with two different treat- 
ments of the absorbing column: with Nn,int = and 
with Nn,int (assumed to be at the redshift of the source) 
allowed to vary freely. We took values of Nn,Gal from 
Stark et al. (1992), although it is probable that the ac- 
tual Galactic absorbing column for each object might be 
slightly different due to the low resolution of the Stark 
et al. survey. In all cases, the opacity associated with 
the intervening column density was modeled by a stan- 
dard Morrison & McCammonn (1983) absorption model. 
The motivation to try both the power law model, as well 
as the less-used logarithmic parabola model (described 
in Giommi et al. 2002 and references therein), was that 
both simple power-law and continuously curving spec- 
tral shapes can reasonably be expected depending on 
the model adopted for synchrotron aging and acceler- 
ation (Leahy 1991, Massaro 2002). Our spectral fits for 
the single power-law model are summarized in Table 3. 
In Figure 1, we show as examples spectra of two objects 
fit with the single power-law model. 

As can be seen, the spectral indices we found range 
from r = 1.7 — 2.9 (corresponding to energy spectral in- 
dices a = 0.7 — 1.9), with twelve of seventeen being in the 
range V = 2.1 — 2.7. These spectral indices are similar 
to the findings of past observations of HBLs by ROSAT 
(Perlman al. 1996a, Padovani et al. 1997), ASCA (Kubo 
et al. 1998) and BeppoSAX (Wolter et al. 1998, Beck- 
mann et al. 2001). Three spectra were found to be flat 
(r < 2). These objects are likely to have synchrotron 
peaks at energies higher than ~ 5 kcV (see also §4). 

In the three objects where repeated observations were 
made, we see significant variability both in flux and spec- 
tral shape (Table 3). In all of those objects, the high- 
est flux state observed also has the hardest spectrum, 
in agreement with the known spectral variability prop- 
erties of BL Lac objects (e.g., Ulrich et al. 1997 and 
references therein). However, the correlation between a 
harder spectrum and higher X-ray flux is not one-to-one, 
as indicated by the three observations of 1ES1959+650. 

3.2. X-ray Spectral Curvature 

In fourteen of seventeen observations, a better fit (sig- 
nificant at the > 99% level according to F-test results; 
see Table 3) was obtained by allowing for the possibility 
of spectral curvature (over and above the default, power- 
law plus Galactic Nh model). This curvature can either 
be intrinsic, or the result of additional absorption due to 
material either within the BL Lac's host galaxy or our 
own. As indicated above, we first investigated the possi- 
bility that additional absorption was present; the result 
of these procedures is given in Table 3. As can be seen, 
in 13 of 14 observations the column required is ^ 20% 
of the Galactic figure. We also investigated the alternate 
possibility of intrinsic curvature, using the logarithmic 
parabola model as well as an alternate model, detailed 
later, where we simply fitted power-laws independently 
in four sub-bands, and obtained results of similar signifi- 
cance. The results of this procedure are given in Table 4 
and discussed below. Importantly, it is impossible to dis- 
criminate a priori between the two models due to their 
similar mathematical forms (see above). It is, however, 
possible to use other information and perform other tests. 

We first inspected the implied columns to check the 



reasonability of the idea that the spectral curvature was 
the result of additional absorbing material. We believe 
that additional absorbing material is less likely for three 
reasons. First, the typical BL Lac's host galaxy is a 
bright elliptical, where columns this high would not nor- 
mally be expected except in extreme cases (Goudfrooij 
et al. 1994). Also, in 1ES1959+650, one of the three 
objects where repeated observations were done, the im- 
plied Nh appears to vary between epochs. This requires 
the absorbing material to have been in a region smaller 
than ~ 1 light-month in size. This would imply densities 
n ^ 10 3 — 10 4 cm~ 3 , which is outwardly not unreason- 
able for the inner regions of the AGN. But the implied 
density would have to climb with decreasing timescale, 
as t~ 3 / 2 , meaning that unreasonable column densities 
would be reached if Nh appeared to vary on timescales 
only 1-2 orders of magnitude smaller, which are not un- 
common for flux variability among BL Lacs (Ulrich et al. 
1997). Finally, there is no indication from optical or UV 
spectroscopy that BL Lac objects would have any signif- 
icant neutral material in the line of sight to the nucleus 
(e.g., Perlman et al. 1996; Rector et al. 2000; Rector 
& Stocke 2001; Kinney et al. 1991; Lanzetta, Turnshck 
& Sandoval 1993; Penton & Shull 1996). This material, 
in principle, can be partially ionized, providing no opac- 
ity in the optical and UV, but absorbing via individual 
edges in the soft X-ray band. This can be ruled out as a 
general property on the basis of an absence of any indi- 
vidual spectral features in the RGS data (§3.3), although 
we cannot exclude it in the spectra of fainter objects. 

This is persuasive but not conclusive evidence that 
these spectra have intrinsic X-ray spectral curvature. To 
test this idea further, we split the XMM-Newton data 
into four sub-bands: 0.5-1 keV, 1-2 keV, 2-4 keV and 4- 
10 keV, and fit each sub-band individually with a power- 
law model. For this exercise we assumed an Nh value 
fixed at Galactic, and used the MOS data only below 
2 keV and both the PN and MOS data at higher en- 
ergies. This was done because of the known issues in 
cross-calibrating PN data with those from the MOS at 
lower energies, especially near 1 keV (see also the first 
paragraph of §3, above). The results of this procedure 
are given in Table 4, and four examples of these fits are 
shown in Figure 3. We note that the overall % 2 values 
for these fits are comparable to those achieved by varying 
Nh; however it is impossible to use an F-test to compare 
them because of the differing ways in which data was se- 
lected for this procedure. As can be seen, the typical 
X-ray curvature seen in our sample is nearly constant, 
even up to the highest energy band. We have char- 
acterized this curvature by looking at the dT/d(\ogE) 
for each object. This information is included in Ta- 
ble 4. As can be seen, our data are consistent with a 
(dT/d(\ogE)) s=a 0.4±0.15 in fourteen of seventeen obser- 
vations (the same ones for which spectral curvature was 
indicated as mentioned above). This is not consistent 
with the observed columns in these objects, which range 
from 3 — 20 x 10 20 cm~ 2 . At such moderate columns, it 
is typical for the spectrum above about 1 keV to be rel- 
atively insensitive to the absorbing column. This is yet 
another argument (albeit indirect) against the presence 
of additional absorbing material explaining the observed 
spectral curvature in these objects. 
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As already noted, we also attempted to model intrin- 
sic curvature by fitting a logarithmic parabola form. The 
information for these fits is also given in Table 4. The 
logarithmic parabola is characterized by both a spectral 
index and a curvature parameter, (3. The spectral form 
is somewhat different mathematically from either of the 
other forms; however, it is essentially similar to the ap- 
plication of a constant dT/d(logE), with the added ad- 
vantage of allowing a consistent fit in the entire band. In 
theory, the value of the curvature parameter /3 should be 
equivalent to the value of (dT/d(logE)) found using sev- 
eral sub-bands. Inspection of Table 4 shows that these 
are indeed similar, albeit with small variations in each 
case. The mean curvature required for the objects where 
this is indicated by the F-tests is essentially identical, in 
fact: (/?) = 0.37 versus (d/T/d(\ogE)) = 0.38. 

Giommi et al. (2002) also found good evidence for 
intrinsic X-ray spectral curvature (which they modeled 
by the same logarithmic parabola form) in about half 
of the HBLs observed by BeppoSAX. The value of {(3) 
they found is consistent with our findings. However, the 
proportion of curved spectra in Giommi et al. (2002) is 
considerably lower than we see. The difference cannot 
be explained by our small sample, as twelve of thirteen 
objects observed by BeppoSAX were also observed by us 
with XMM-Newton, and of those, for six the best repre- 
sentation of their BeppoSAX spectrum was a logarithmic 
parabola plus fixed absorption (although n.b., Giommi et 
al. 2002 did not fit a variable absorption model for any 
object). We believe it is much more likely that the reason 
for the larger fraction found here lies in the higher sig- 
nal to noise and resolution of these spectra compared to 
those from BeppoSAX. This is supported by the fact that 
if one examines the variable Nu fits, the 6 objects where 
Giommi et al. (2002) required a logarithmic parabola 
are best fitted by significantly higher absorbing columns 
than those which Giommi et al. (2002) was able to fit 
adequately with a single power law. 

Our conclusion from the above is that we believe the 
most appropriate model for the spectral curvature we 
see is intrinsic, rather than due to additional absorption. 
Given the similar values of x 2 , our primary motivation 
for concluding this is the indirect lines of evidence pre- 
sented above regarding the implausibility of additional 
absorbing columns under these circumstances. We will 
proceed on this basis in the upcoming discussion. 

To close this discussion of intrinsically curved spectra, 
it is useful to compare the fits and shape of the spec- 
trum achieved for a given object under a simple power 
law model (with Galactic Nh) with those implied by the 
logarithmic power law form. We show this comparison in 
Figure 4, for one of the spectra of 1ES1959+650, specif- 
ically OBSID 0094383301. There comparison is shown 
in two forms: at top, we show the comparison in terms 
of the detected spectrum (folded into the instrumental 
response), while at bottom it is shown in terms of emit- 
ted energy. The plots at top enable us to see the real 
differences in the quality of the fit: the simple power 
law model (Figure 4, top left) clearly overprcdicts the 
observed spectrum at E > 1 keV, and slightly under- 
predicts the lower-noise MOS data at 1-1.5 keV. But at 
higher energies it is the PN data which have the greater 
ability to distinguish between the two models, and the 
simple power law model overpredicts the data at energies 



higher than about 5 keV. As can be seen, the logarith- 
mic power law model (Figure 4, top right) fits the data 
much better. The plots at bottom enable us to see the 
form of the curvature: note that above ~ 1 keV essen- 
tially no curvature would be seen, if the spectrum fol- 
lowed a simple power law model (Figure 4, bottom left). 
However, under the logarithmic parabola model (Figure 
4, bottom right), we see that there is in fact significant 
curvature seen at higher energies, albeit of a mild, grad- 
ual form. This is presumably why this curvature was 
not seen by earlier satellites, which suffered from cither 
a much smaller spectral range (for all satellites except 
BeppoSAX) and/or much lower signal to noise (all other 
satellites). 

3.3. Spectral Line Features in BL Lacs? 

We inspected each spectrum thoroughly for spectral 
features. We did not find convincing evidence for line 
features in any of these spectra. To quantify this, both 
the EPIC and RGS data for the 12 highest signal-to- 
noise spectra were searched for absorption features in 
XSPEC (we did not search datasets where either instru- 
ments were turned off and/or the statistics were inad- 
equate). The best- fit broadband model (Table 4) was 
used as the seed model in each case, and the RGS data 
were grouped in bins with a minimum of 20 photons 
apiece (following our practice for the EPIC data; §3.1). 
We then artificially set a line feature of 5 eV physical 
width, and covering fraction of 0.5, and then stepped it 
through the range 0.5 to 4.5 keV in increments of 0.01 
keV using STEPPAR, and assessed the significance of 
features using the Ax 2 (A 3cr feature has Ax 2 = -9.21, 
while a 4cr feature has A^ 2 = —18.42). Higher ener- 
gies were not searched due to the small effective area 
and poor statistics. Once features were found, we then 
fixed the energy and allowed XSPEC to converge on a 
best-fit width. This procedure found six narrow fea- 
tures, as shown in Table 5, in four of the datasets, but 
in most (1ES 0120+340 OBSID 0094382101, 1ES0323 
OBSID 0094382501, 1ES 0414+009 OBSID 0094383101, 
1ES0647+250 OBSID 0094380901, 1ES1133+704 OB- 
SID0094170101, 1ES1255+244 OBSID 0094383001, 1ES 
1553+113 OBSID 0094380801, 1ES1959+650 OBSID 
0094383501), no features were found. 

Are these features real? As can be seen, all but one are 
between 3-4 a, with one barely over 4a. For Gaussian 
statistics, P(3a) = 0.0027, while P(4cr) = 0.000318. As 
can be seen from Table 3, the typical EPIC dataset in our 
sample is ~ 300 channels, as is the typical RGS dataset 
(n.b., the number of channels is not significantly greater 
due to the relatively poor statistics resulting from the 
short exposure times). Thus the result of 5 features at 
3c and 1 feature at > 4cr in ~ 3600 spectral channels 
is consistent with what one would expect from statistical 
noise (respectively 10 + 3 features at > 3cr and 1 + 1 
feature at > 4a). Moreover, as noted in Table 5, none 
of the features found in any of the spectra are present 
in the spectra from all the EPIC and RGS instruments. 
Thus, we believe all are the result of random noise and 
therefore not real spectral lines. The limits shown in 
Table 5 can be taken as representative of the detection 
limits of our data, which range from ~ 5 up to ~ 50 eV 
depending on count rate. 

This conclusion differs from findings based on data 
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from earlier instruments. We are confident of our result 
based on the relatively high signal-to-noise of the XMM- 
Newton spectra. We note, however, that narrow, small 
equivalent width absorption lines such as those inferred 
(at intermediate redshifts) from Chandra grating obser- 
vations of BL Lac objects by Nicastro et al. (2002; see 
also Cagnoni et al. 2004), or those due to ionized oxygen 
in the hot halo of our Galaxy (e.g., Nicastro et al. 2002, 
McKernan, Yaqoob & Reynolds 2004), are not excluded 
by the data in-hand. The determination of presence or 
absence of such features require much higher S/N data 
than are available in our sample. 

It is worth noting that a similar absence of strong 
spectral features was also found in another, independent 
study of the XMM-Newton X-ray spectra of four of the 
five BL Lacs where BBXRT and ASCA spectra appeared 
to show these line features, specifically H1219+301, Mkn 
501, H1426+428 and H0548-322 (Blustin et al. 2004). 
Those authors reported that if any of the previously 
found absorption features were indeed real, they had 
to represent a transient phenomenon. They found this 
particular conclusion unlikely, ruling it out at 93% con- 
fidence based on the assumption that the (multiple) 
XMM-Newton observations of those objects represented 
random time instances. The objects in the Blustin et al. 
(2004) sample are not discussed here, and moreover, our 
integration times are smaller so we have used primarily 
the EPIC data and not the RGS. However, our results 
arc in the same vein. With the addition of our data, it 
appears unlikely that the absorption features found with 
ASCA and BBXRT were transient phenomena. 

The most likely explanation of the previously reported 
line features is a combination of mis-calibration of the 
previous instruments and the use of overly simplistic 
spectral models. As an example, a deficit of low energy 
counts in the ASCA data has been previously noted in 
reference to the data for NGC 5548 (Iwasawa, Fabian, 
and Nandra 1999); the simultaneous ASCA and ROSAT 
PSPC observation implied significantly higher flux be- 
low 0.5 keV than any reasonable spectral extrapolation 
of the ASCA SIS data would allow. Another reason is 
that a spectral curvature of exactly the type discussed in 
this paper would, at low resolution and lower signal-to- 
noise, mimic an absorption feature when the underlying 
spectrum was assumed to be a simple power law. How- 
ever, the confirmation that an absorption feature was 
measured would require data below the energy of such a 
putative feature, and such data were not always available 
(as was the case for the Einstein SSS as well as SSS + 
MPC data). We therefore suggest that the claims made 
by earlier workers were a byproduct of the particular 
simple power law model used, while the effect was due 
to gradual spectral curvature and/or instrumental effects 
rather than absorption features from material along the 
line of sight. 

4. MODELING THE X-RAY AND BROADBAND 
SPECTRAL CHARACTERISTICS 

We have constructed broadband spectral energy dis- 
tributions (SEDs) for each object using the EPIC data, 
in conjunction with the data from the OM and the UM- 
RAO flux nearest in time to the XMM observation. For 
objects which were too faint in either the optical or radio 
at the time of observation to register a positive detection 



on either the OM or UMRAO, we have used POSS plates 
to obtain an average optical flux and radio fluxes from 
the literature (Perlman et al. 1996b and sources therein) 
for the SED. We fit to these data a simple parabolic 
model, in order to obtain a peak frequency, f pea k- We 
did not assume a standard synchrotron spectrum for this 
procedure because the small amount of data we have do 
not adequately constrain these models (see e.g., Leahy 
1991, Pacholczyk 1970). For a few objects, where the 
initial fit was not good, we had to modify a single non- 
simultaneous point to find a convincing value of v pea k . 

In Table 6 we summarize the broadband spectral data 
and v pea k for every object in our sample. As can be 
seen, our objects range between log v pea k ~ 14.5 — 18, as 
expected for HBL-type and intermediate objects. A few 
objects, noted in Table 6, have very flat X-ray spectral 
index and/or significant evidence of v pea k > 10 18 Hz. 
Our procedure yielded is pe ak values considerably in excess 
of 10 18 Hz for all these objects (see also §3.1), but in 
Table 6 we list them as having v pea k — 10 18 Hz because 
our data are not able to test for a higher peak frequency 
given their frequency space coverage. 

We then proceeded to test for correlations be- 
tween broadband spectral properties and both T and 
dT/d(logE). We then tested each for correlations using 
standard statistical tests (Spearman's rank correlation, 
p, and Pearson's correlation coefficient, r). These plots 
and their physical implications are the main subject of 
this section. 

In Figure 5 we show plots of T versus a ro (left), 
a ox (middle) and v pea k- As can be seen, our data 
show significant anti-correlations between T and v pe ak 
[r = -0.741, p = -0.781, P = 0.18%) and between P 
and a ox (r = 0.609, p = 0.634, P = 1.1%) but no signifi- 
cant correlation between T and and a ro (r = 0.036, p — 
— 0.05, P = 84.0%). This constellation of results is con- 
sistent with the a x - v pea k relation of Padovani et al. 
(1997), when one considers that the range of v pea k val- 
ues we cover is only ~ 3 decades. 

4.1. Statistical Particle Acceleration? 

The observed spectra exhibit curvature over a fre- 
quency range from 0.5 keV to 10 keV, in the sense 
of a gradual steepening with increasing energy. As- 
suming that the observed emission is synchrotron ra- 
diation, this implies that the electron distribution re- 
sponsible for this emission is not the typical power law 
predicted by particle acceleration schemes (for a recent 
review see Gallant 2002). Instead, the electron spec- 
trum has to gradually curve downwards in a way that 
produces the observed spectra. Given that X-ray syn- 
chrotron emission requires in situ particle acceleration 
due to the short radiative lifetimes of the emitting par- 
ticles (T sync = 1.2 x 10 3 B G 3/2 E^ 2 S-^ 2 s; Rybicki & 
Lightman 1979), the ultimate physics behind curved X- 
ray spectra must be connected on a very basic level with 
the nature of particle acceleration in these systems. Here 
and in the next subsection we discuss possible solutions 
to this problem. 

In Massaro et al. (2004b) a log-parabolic spectral cur- 
vature was fit to the X-ray spectrum and broadband SED 
of Mkn 421 in several flux states. Those authors mod- 
eled the log-parabolic curvature via a statistical particle 
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acceleration process under which a particle i with energy 
7i in a given region has a finite probability (< 1) of being 
accelerated, where the probability is given by pi = g/lf, 
where g, q are positive constants. They show in their 
paper that this leads analytically to a log-parabolic cur- 
vature form (see their §6 and equations (10)-(18)), and 
predicts a linear correlation between the X-ray spec- 
tral index, r x , and dT/d(logE). In Figure 6, we show 
plots of Vpeak versus dT / d(logE) (top) as well as T x ver- 
sus dT /d(logE) (bottom). As can be seen, neither of 
these plots shows a significant correlation (v pea k versus 
dT/d{logE) : r = 0.615, p = 0.393, P = 11.6% ; T x ver- 
sus dT/d{logE) : r = -0.24,/) = 0.014, P = 95.3%). 
The results are similar for tests of v pe ak and T versus /?. 
Thus our data do not support the model of Massaro et 
al. (2004b). One possible reason why this is so is that 
Massaro et al. (2004b) model particle acceleration only, 
and do not include losses. Moreover, the statistical par- 
ticle acceleration is still assumed to be time-invariant at 
all locations in the jet and hence the issue of variability 
is not addressed. 

4.2. Spectral Curvature: a Signature of Episodic 
Particle Acceleration? 

The results of §4.1 prompt us to look for another expla- 
nation that takes into account the known observational 
properties of jets and particularly blazars. Another mo- 
tivation for this is the recent result of Perlman & Wilson 
(2005) that standard, continuous injection models of par- 
ticle acceleration consistently overpredict the X-ray flux 
of the M87 jet by large factors (up to ~ 100 at 1 keV, 
varying with energy roughly as E 0A ). The interpreta- 
tion advanced in that paper was one of a position- and 
energy-variable (but not time-variable) filling factor for 
particle acceleration. The Perlman & Wilson work, how- 
ever, does not explore the issue of spectral curvature, 
as the statistics in most components of the M87 jet are 
inadequate for this purpose; nor does it explore issues re- 
lated to variability, which are known to be of paramount 
importance for BL Lac objects. 

In this subsection we address these issues directly, by 
demonstrating that curved electron particle spectra can 
be produced if the particle acceleration is episodic. Im- 
portantly, this model assumes only a time-variable par- 
ticle acceleration which in an integrated sense should be 
indistinguishable from a a spatial average. Unlike the 
case of M87, here we have no spatial information regard- 
ing the distribution of jet X-ray emission and/or particle 
acceleration. Instead, we know that variability is an im- 
portant characteristic of blazar jets and thus we optimize 
our model for this case. We explore these subjects more 
deeply in a later paper (Georganopoulos & Perlman, in 
preparation) that fully outlines a model considered in 
brief form below. 

Consider a zone, possibly a shock, where electrons with 
Lorentz factor 70 are injected at a rate of Q electrons 
per second and accelerated to higher energies. Following 
Kirk, Rieger, & Mastichiadis (1998), the kinetic equation 
describing particle acceleration is 



dn 



d_ 



1 



■0Y 



+ 



QS(7-7o), (1) 



where n(7, t) is the electron energy distribution, l/t a 



is the particle acceleration rate, \jt esc is the escape rate 
of particles from the acceleration region, and 
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with o~t the Thomson cross-section and U the total en- 
ergy density of the magnetic field and ambient photons. 
Assuming that the injection started at time t = 0, the 
electron energy distribution after time t will be a power 
law with electron index s = 1 + t acc /t esc up to a Lorentz 
factor 
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where jmax — 1/ ' Pt aC c is the maximum energy electrons 
reach asymptotically if the particle acceleration mecha- 
nism operates for t > ln[(7 TOax - 7o)/7o]£a CC - The time 
t required for electrons to be accelerated up to Lorentz 
factor 71 is 
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If we assume that particle acceleration is characterized 
by a typical timescale T, then the electron distribution 
at any time < t < T is a power law that cuts off at a 
Lorentz factor l(t), reaching its maximum Lorentz factor 
7i(T) at t = T. Particle acceleration ceases to operate 
at t = T, and the total time r over which electrons of 
Lorentz factor 7 are provided by the particle acceleration 
mechanism is 
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The time r is also the time during which an electron of 
Lorentz factor 7 radiates, as long as the radiative cooling 
timescale t rad = l//?7 = t acc -f max /-y satisfies the condi- 
tion t rad <C t. In other words, r is the total time the 
accelerated electron distribution is non-zero at electron 
Lorentz factor 7; note that, as is apparent from Eq. (5), 
the condition t rac i -C r can always be satisfied for 71 
sufficiently close to l max - If the light crossing time of 
the emission region or the integration time of our obser- 
vations is greater than T, then we effectively observe a 
time averaged electron distribution (77,(7)) which is the 
product of a power law multiplied by a logarithmic term 
representing the time interval this power law is available 
at each energy: 
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Using the the 5- function approximation (i.e., an electron 
of Lorentz factor 7 produces synchrotron photons only 
at the critical synchrotron energy) to calculate the syn- 
chrotron emission, we obtain an observed photon spec- 
trum 
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where E\ and E max are the synchrotron energies pro- 
duced by electrons of Lorentz factor 71 and jmax re- 
spectively. In Figure 7 we show the photon spectrum 
and photon index as a function of energy, for the fol- 
lowing parameter choices: s = 3.0, E max = 100 keV and 
Ei = 50keV. As can be seen, the spectrum steepens grad- 
ually, with r increasing at a fairly constant rate up to ~ 
few keV, and then dramatically at E > 10 KcV. If our 
interpretation of the steepening is correct, we expect to 
see an increasing curvature in future observations above 
the 10 keV upper limit of XMM-Newton, as predicted 
in equation (7) (see Figure 7). One might also expect 
that in other objects, e.g., LBL, this cutoff might occur 
at lower energies. In LBL, however, we do not observe 
such steep values of T, but rather, flat X-ray spectra are 
seen (r < 2; see Padovani et al. 2001, 2004), which 
are typically interpreted as being dominated by inverse- 
Comtpon emission (Padovani et al. 1997, Lamer et al. 
1996). In practice, it may be difficult to observe this 
steepening because as the curvature of the synchrotron 
component increases, the spectrum is essentially cutting 
off and at some point the observed spectrum will begin 
to be dominated by the onset of the Compton compo- 
nent. What one might expect to see in practice, then, 
would be concave spectra (i.e., increasing T) up until the 
two components cross, with an inflection point at their 
crossing point and then a constant value of T (< 2) at 
higher energies. One final thing that should be noted 
here is that above 10 keV, almost all the HBL spectra 
that currently exist are for objects that are currently in 
a flaring state, where we would not expect to see signif- 
icant steepening because the particle acceleration is still 
at or near its maximum. Also, there arc no currently 
available instruments that would be sufficiently sensitive 
at E > 10 keV to search for further steepening in the 
spectrum of HBLs. 

The two energies E\ and E max appearing in Equation 
(7) can be derived through spectral fitting. However, 
their values alone are offering little insight to the phys- 
ical conditions in the source, because they depend on 
a series of parameters (B, t acc , 70, T). As we discuss 
in §5, additional information from frequency-dependent 
variability can be used to constrain the physical state of 
the source. 

5. CONCLUSIONS 

We have presented 17 high-signal-to-noise XMM- 
Newton X-ray spectra for 13 BL Lac objects. These 
spectra do not show any spectral lines, contrary to ex- 
pectations from some previous observations. We believe 
the reason for this discrepancy is due to a combination 
of overly simplistic spectral models and/or insufficiently 
precise calibration of previous instruments. As a result 
the intrinsic X-ray spectra of BL Lac objects appear to 
be quite featureless, giving no sign as to the physics of 
any region within the central engine, other than the jet 
itself. This is quite different from all other non-blazar 
type sources, where some line features, in either emission 
or absorption, are seen. This can be interpreted in the 
context of unified scheme models (e.g., Urry & Padovani 
1995), under which the jet emission is Doppler boosted 
because it is seen at small angles to our line of sight. It 
is, however, difficult to address the issue of the necessary 
viewing angles to produce the degree of fcaturelessness 



seen in these spectra, particularly given that at present 
our knowledge concerning the X-ray spectra of the par- 
ent population of FR 1 radio galaxies is limited. Fu- 
ture work will provide much information on this subject, 
but currently the best information available (based on 
luminosity function estimates) suggests viewing angles 
6 ~ 10 — 30 and Lorentz factors r ~ 5. This would yield 
Doppler boosts to the jet's apparent luminosity of hun- 
dreds to thousands, as detailed in e.g., Urry & Padovani 
(1995). 

Our data are best fit by curved spectra which grow 
steeper in an approximately logarithmic fashion across 
the XMM-Newton band. From an observational point 
of view, a logarithmic curvature is not a new concept, 
having been first introduced by Landau et al. (1986) 
for radio-optical observations of blazars, and then re- 
vived for X-ray spectra by Giommi et al. (2002). Our 
analysis is, however, the first to demonstrate that in- 
trinsic curvature is to be preferred over the alternate 
model of additional absorbing material along the line of 
sight. From the theoretical point of view, logarithmically 
curved spectra require explanation, as most of the stan- 
dard models of synchrotron emission produce spectra of 
a constant power-law slope (see, e.g., Leahy 1991 for a 
review). Massaro et al. (2004b) suggested statistical ac- 
celeration, whereby a given particle in the distribution 
has a probability of being accelerated, as one possible 
model for this type of spectral shape. The predictions of 
that model are not borne out by our analysis, however, 
and a statistical model also does not take into account 
the variable nature of BL Lacs. It was with this moti- 
vation that we considered an episodic model for particle 
acceleration, which explicitly considers a modification to 
the kinetic equation and thus includes both losses and a 
model of variability. We showed that this model fits our 
data reasonably well and predicts continued steepening 
at higher energies, until the onset of the IC component. 

Episodic particle acceleration may be related to the 
broadband variability observed in TeV sources like 
1959+650 (Krawczynski et al. 2004), where both the 
power and peak frequency of the synchrotron spectrum in 
X-rays exhibit an increase by more than a factor of ~ 10, 
while the optical flux remains practically constant. This 
can be easily accommodated in our scheme through an 
increase of the characteristic time T that particle acceler- 
ation operates. In this case 71 and E\ will both increase, 
resulting in a more powerful X-ray synchrotron spectrum 
that peaks at higher energies. The situation will be quite 
different at lower frequencies; the curvature of the spec- 
trum produced by the logarithmic term will be present 
down to photon energies for which the radiative cooling 
times t r ad becomes larger than the characteristic parti- 
cle acceleration time T. Below this energy the electron 
energy distribution and the corresponding synchrotron 
spectrum will be a simple power law. At these lower en- 
ergies an increase in T can only affect the amplitude of 
the observed spectrum and even that only to the extent 
that such an increase results in an increased fraction of 
time that particle acceleration is operating (often called 
the duty cycle). If, for example, T describes both the 
time acceleration is on and off, then as T increases the 
average number of particles injected remains constant, 
resulting in a constant flux at lower frequencies, while 
the spectrum changes dramatically at high frequencies 
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due to the increase of both 71 and E\ . It is very encour- 
aging that our observations of 1ES1959+650, one to three 
months after flare observed by Krawczynski et al. (2004), 
show significant spectral curvature. This indicates that 
the characteristic timescales on which episodic acceler- 
ation might operate is ^ 1 month, which is consistent 
with the findings of long-look campaigns (Perlman et al. 
1999, Tanihata et al. 2001). Such a timescale is also con- 
sistent with Doppler-boosted versions of the flare found 
recently in the jet of M87 (Harris et al. 2003, Perlman 
et al. 2003), requiring only values 5 ~ 10 (as opposed to 
the more modest 6=1 — 2 required for M87). Indeed, 
the fact that we see logarithmic curvature in 14 of 17 of 
these spectra indicates that the typical ratio of t/T ^ 1; 
i.e., the amount of time in which flares are seen is not 
appreciably smaller than the amount of time in which 
more gentle or no variability is seen. This is not inconsis- 



tent with the data from the ASCA long-look campaigns 
(Tanihata et al. 2001). 
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Fig. 1— Two examples of the XMM EPIC spectra. At top, Mkn 180 (OBSID 0094170101) and at bottom, 1ES1959+650 (OBSID 
0094383301). In both panels, the EPN data are the top line and the data from the MOS detectors are below. The model shown represents 
the power law plus variable Njj fit. See Sections 2 and 3 for discussion. 
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Fig. 2. — Radio light curves for 1ES 1553+113 (top) and 1ES 1028+511 (bottom), obtained at the University of Michigan Radio 
Astronomy Observatory. Data at 14.5, 8.0 and 4.8 GHz are denoted by crosses, circles and triangles respectively. Daily averages of the 
data are shown, with the bars denoting la error estimates. As seen, both sources vary considerably during the radio observations. Sec 
Section 2.1-2.2 for discussion. 
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Fig. 3. — Four examples of fitting the spectra of BL Lacs in our sample to smaller spectral ranges. All four of these objects show steeper 
spectra at higher energies, with dT/d(logE) Ri 0.4. See Section 3.3 for discussion. 
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Fig. 4.— Comparison of model fits for the XMM spectrum of 1ES1959+650 (OBSID 0094383301). In the panel at left, the data were 
fitted with a power-law plus Galactic N(H) model, while at right, the data were fitted with the log-parabola plus Galactic N(H) model. 
Note that the fit in the left panel is adequate between 1.5-5 keV, but cannot account for the curvature seen both at energies < 1.5 keV and 
> 5 keV. At top, the comparison is shown in terms of the detected spectrum, with each plot showing received spectra plus folded model 
for the MOS 1 & 2 (bottom lines) and PN (top lines), as well as residual. At bottom, we show the unfolded energy spectrum. These plots 
show that curvature is present in this spectrum up to the highest energies seen by XMM. 
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Fig. 5. — Plots of the X-ray spectral slope T, versus broadband spectral index a ro (top left) and a ro (top right), as well as the SED peak 
frequency u pea ^ (bottom left). For easy cross-referencing with the Tables, we note in these plots the identity of each point. See Section 4.1 
for discussion. 
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Fig. 6.— Plots of dT/d{logE) versus T (top) and /^p ea fc (bottom). For easy cross-referencing with the Tables, we note in these plots the 
identity of each point. See Section 4.1 for discussion. 
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Fig. 7. — Upper panel: The photon spectrum (solid line) resulting from an episodic injection of electrons in the particle acceleration 
process. Here s = 3.0, E max = 100 KeV, and E\ = 50 KeV. Note that E\ is the maximum emitted energy the source is reaching in this 
episodic mode of injection, while E max is the emitted energy the source would reach under steady injection conditions. For comparison, 
a power law emission corresponding to an electron distribution with s = 3.0 is also shown (broken line). Lower panel: the photon index 
T that corresponds to the emitted photon spectrum. Note the gradual steepening of the spectrum, manifested by an increase of T with 
energy. The dotted lines represent the energy range 0.5 — 10 keV covered by XMM. See Section 4.2 for discussion. 
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39.5 




0094373501 


2003-02-09 


3700 


Medium 


PartialW2 


12.5 (3.80) 3 


2000 


Medium 


Small 


43.3 



x This observation was affected by proton flaring 

2 Observation was attempted again in the same orbit 

3 Count Rates are raw count rates in the effective exposure interval, while values in parentheses are rates after pile-up correction. 
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TABLE 2 
Supporting observations 



1ES Name OBSID OM Band F(UV) -F max (rad) 









(10 10 erg cm z s 1 A i ) 


(Jy) 


0033+595 


0094381301 


UVW2x2 


< 3.3 1 


0.3 


0120+340 


0094382101 


UVW2x5 


8.9 ± 0.9 


0.1 


0145+138 


0094383401 


U VW2 x 4 


< 3.3 1 


0.2 


0323+022 


0094382501 


UVW2x5 


14.4 ± i.o 


0.3 


f)347— 1 21 


0094381 1 01 


VJ V VYiAl 


7.6 ± 2.0 


2 


041 zi4-nno 


ooqa^&^i ni 


T TVW9 v \ 

U V VV Zr A O 


io.o in u. ± 


1 

u. 1 


nfi/i i i oc;n 
U04 / -f-ZOU 


nno/i Qcnoni 
UUy4ooUyUl 


U V W Z X 1 


< 0.0 


n 1 
U. 1 




0094382901 


UVW2 x 5 


< 3.3 1 


0.1 


1028+511 


0094381801 


U VW2 x 5 


18.2 ± 1.6 


0.2 




0094382701 


UVW2 x 4 


21.9 ± 0.7 


0.2 


1101-232 


0094380601 


U VW2 x 4 


10 ± 3 


0.3 


1133+704 


0094170301 


N/A 


N/A 


2 




0094170101 


Ux4 


272 ±8 


N/A 3 


1255+244 


0094383001 


U VW2 x 5 


3.2 + 1.5 


0.2 




0094383201 


U VW2 x 4 


< 3.3 1 


0.2 


1553+113 


0094380801 


U VW2 x 5 


185 + 5 


0.8 


1959+650 


0094380201 


U VW2 x 5 


35.0 + 2.0 


0.6 




0094383301 


U VW2 x 5 


25.3 ± 2.0 


0.6 




0094373501 


U VW2 x 5 


28.4 + 2.0 


0.6 



1 Object not detected in the OM observation; value is a 2cr upper limit 

2 Below detection threshold in all UMRAO observations 

3 No UMRAO observations made as this object was originally part of another program 



TABLE 3 
BL Lac Spectral Fit Parameters 



1ES Name 


Obsld 




Fit parameters, N(H) Fixed 






Fit parameters, N(H) Free 








N H ,Gal 


T F2-10keV 


X 2 /dof. (dof.) 


N H ,tot 


r F 2 _ 10kcV X 2 /dof. (dof.) 


F-tcst 






10 20 cm" 2 


erg cm~ 2 s _1 




10 20 cm" 2 


erg cm 2 s 1 


Probability 



00^+^9^ 


0094^81 ^01 F 


42.4 


9 9+0.02 
z ' z -0.03 


1.2 




10" 


11 


1 4.9(9%%} 


f?n q+3.8 


9 47+0 0^ 


1.1 x 10" 


n 


1 01 (9R7\ 


<? in -10 


01 20+340 


0094382101 


5.14 


2 16+ ' 01 


9.7 


x 


10" 


12 


1 70('527 > ) 


9.9±0.6 


2 34+0- 02 
z -° -0.03 


8.9x10" 


12 


1 32('526 > ) 


< io- 10 


01 4^-Ul ^8 

U14iJT lOO 


00QzL383zL01 


5.10 


9 07 + 19 


2.5 




10" 


13 




99 1+8.0 


9 no+0.38 


1.6x10" 


13 


1.24(87) 


3.2 x 10 - 4 


0393-1-099 


ooqa3&9^oi 


R 7 A 
o. ( 4 


9 c-4+0.02 


9 




"1 o - 


12 


1 0^91 3^ 


11 9±1 2 


9 fi7+0 OR 


1 9 x 10" 


12 


n 0^91 o\ 


1 3 x 10" ^ 


0347—121 


0094381 101 


3.64 


i 77 +0.02 
1 - ' ' -0.01 


1.8 


x 


10" 


11 


1 ISfSSO 1 ) 


5.0±0.6 


1 82+0 03 


1.8x10" 


11 


1 lfi('52Q s ) 


2.6 x 10 -3 


0414+009 


0094383101 


10.3 


9 47 +0.02 

z - 4 '-o.oi 


9.3 


X 


io- 


12 


1.20(499) 


13.1±0.6 


2.57±0.03 


8.9x10- 


12 


1.09(498) 


< IO" 10 


0o4l +250 


0094380901 


12.8 


9 ^9 _l_ n n9 
z.bz ± 0.02 


5.0 


X 


10 


12 


1.30(2 to ; 


19 -°±i'.i 


2.8o±0.05 


5.0x10" 


12 


1.05(2 //) 


^ in — 10 
< 1U 


1028+511 


0094381801 F 


1.16 


-0.03 


1.2 


X 


io- 


11 


1.20(107) 


3.77±1.1 


2.61±0.08 


1.0x10" 


11 


1.04(106) 


< IO" 10 




0094382701 


1.16 


2.28 ±0.01 


1.2 


X 


io- 


11 


1.41(600) 


5.2±0.4 


2.46±0.02 


9.7x10- 


12 


0.98(599) 


< IO" 10 


1101-232 


0094380601 F 


5.76 


9 90+O.O2 

z - zo -o.oi 


2.3 


X 


io- 


11 


1.37(682) 


10.9±0.6 


2 40+ ' 02 
z - 4u -0.03 


2.2x10- 


11 


1.05(681) 


< IO" 10 


1133+704 


0094170101 


1.42 


2.40 ± 0.02 


9.8 


X 


io- 


12 


1.21(361) 


1.8±0.7 


2.42±0.04 


9.7x10- 


12 


1.20(360) 


8.4 xlO- 2 


1255+244 


0094383001 F 


1.26 


1 99+ 01 


5.1 


X 


io- 


12 


1.18(380) 


2.9± 2 2 9 6 


1.83±0.12 


4.8x10- 


12 


1.02(379) 


< IO" 10 




0094383201 


1.26 


1 76+ 07 
x - ,D -0.06 


7.5 


X 


io- 


12 


0.97(134) 


4.6±0.8 


9 -i i +0.04 
z - 11 -0.03 


4.8x10- 


12 


0.97(133) 


N/A 


1553+513 


0094380801 


3.67 


2.38 ±0.01 


3.5 


X 


io- 


11 


1.26(815) 


5.7±0.4 


2 46+ ' 01 
z - 4D -0.02 


3.4x10- 


11 


1.17(814) 


< io- 10 


1959+650 


0094380201 F ' a 


10.1 


1.70 ±0.05 


2.9 


X 


io- 


10 


1.44(58) 


13 3+ 31 


l-80-o.ii 


2.8x10- 


10 


1.41(57) 


0.275 




0094383301 F 


10.1 


2 10+ - 02 


7.7 


X 


io- 


11 


1.51(514) 


17.9±1.0 


9 9Q+0.02 
^•^ a -0.03 


7.2x10- 


11 


1.08(513) 


< io- 10 




0094383501 


10.1 


2.02 ± 0.02 


6.9 


X 


io- 


11 


1.64(863) 


16.9±0.6 


9 1O+0.02 

Als -o.oi 


6.5x10- 


11 


1.13(862) 


< IO" 10 



Note. — MOS data used in 0.5 - 10.0 keV range PN in 1.1 - 10.0 keV range. 

a No PN data available, fit parameters are for MOS data only. 
F Most of this observation is contaminated by proton flaring. 



TABLE 4 

Parameters for BL Lac Spectral Curvature 



4-band T-Fitting 12 Logarithmic Parabola™ 



Lhio iNamc 


UDsld 


J- 0.5-1.0fceV 


1 1.0-2.0fceV 


1 2.0-4.0fceV 


1 4.0-10. OkeV 


XZ/aoi. (aoi.j 


al / (lylogrj ) 


1 


a 
P 


^,0 JA^f ( A^f \ 

XZ/dol. (aoi.j 


0033+595 


0094381301 


i-49ti;° 9 I 


2 03+ ' 18 


2 S4+ ' 12 
Z "J^_0.12 


q q o-hO . 22 

Z " 3J — 0.22 


1.04(187) 


0.58 


1 71 01 

lm ' -0.09 


-o.67i8:l 2 


1.10 (287) 


0120+340 


0094382101 


i.94i8;i 


2 14 +0.08 
z - ±4 -0.08 


2 19+° 1 
z - ±M -0.09 


q on +0.22 
z - oz -0.16 


1.16(383) 


0.4 


2.05 ± 0.02 


-0.30 ±0.04 


1.37 (526) 


0145+138 


0094383401 


0.25l§:» 


9 OO+0.72 


2 ,4+ - 91 
z -3 4 _0.71 


n , c + 1.12 

0.16I 2 76 


1.34(77) 


-0.76 


i.84i8: 2 7 


-1 26+ ' 56 
1 - zu -0.68 


1.27(87) 


0323+022 


0094382501 


2 29+ ' 17 
z - zy -0.17 


2 6 +0 ' 15 
Z - D -0.15 


9 c 9 +0.18 
z - oz -0.19 


9 SQ+ ' 53 
z - 8y -0.38 


1.03(156) 


0.47 


z -48l8:8! 


n 99 +0.09 
u - zz -0.08 


0.96(212) 


0347-121 


0094381101 


J -°-0.13 


1 o 9 +0.09 
1 - oz -0.09 


1 q+01 
1 - a -0.11 


! 7^+0.13 
!• 'O-0.17 


1.15(398) 


0.22 


1 7c -+0.03 
1 - ,o -0.02 


-0 05+ 05 


1.18(529) 


0414+009 


0094383101 


24 +0.11 
z -**_0.11 


2 45+ ' 09 
• 2 -*->_0.09 


9 cri+0.09 
z - J1 — 0.09 


9 S9 +0.18 
• 2 -° z _0.18 


1.15(309) 


0.34 


2.41 ±0.02 


-0.16 ±0.04 


1.11(498) 


0647+250 


0094380901 


2 44+ - 16 
z - 44 -0.16 


2 56 +0 13 
z - do -0.13 


9 7 +0.17 
z - ' -0.18 


2 - 89 -o!38 


0.89(197) 


0.45 


9 tr + 0.03 
z - d -0.04 


-0.39 ±0.08 


1.10(277) 


1028+511 


0094382701 


9 9 o+0.07 
Z - ZJ -0.07 


9 9 o+0.07 
z - ZJ -0.06 


2 4 n+0.09 
Z -* M -0.09 


2.48l°;i 7 


0.93(425) 


0.36 


2.20 ±0.01 


-0.30 ±0.03 


1.00(599) 




0094381801 c 














— 2.42 ±0.04 


-0 19+° 10 

u - la -0.12 


1.13(106) 


1101-232 


0094380601 


2-06±g;ii 


9 1 o+0. 09 
z - lo -0.09 


9 oc+0.08 


9 OQ+0. 16 


0.93(472) 


0.41 


2.09 ± 0.02 


n o 9 +0.03 
u - oz -0.04 


1.05(681) 


1133+704 


0094170101 


9 o+0.14 
z -°-0.14 


9 91 +0.14 
z - zi -0.14 


2 S4+ - 14 
z -" 54 -0.14 


2 47+0.28 
z - 4 '-0.27 


1.28(223) 


0.14 


9 oQ+0.03 
z - c>y -0.02 


-0 01 + 05 


1.21 (360) 


1255+244 


0094383001 


1.88 ±0.13 


1.99 ±0.11 


2 06+ ' 12 
Z - UD -0.13 


2.24lg; 2 l 


0.97(283) 


-0.5 


1 Q0+ 003 
i - 9u -0.02 


-0.22 ±0.05 


1.05(379) 




0094383201 


2 29+ ' 48 
z - zy -0.49 


1.70 ±0.37 


1 71 +0.41 
i -' 1 -0.45 


l.08t°;« 


0.84(112) 


0.35 


1 76+°' 10 
x - ' D -0.10 


-0 01+ 018 
u - ui -0.21 


0.98(133) 


1553+113 


0094380801 


9 oc+0.07 


9 9 rr + 0.07 
z - zo -0.07 


2 42+O.O6 
z -^ z -0.06 




1.11(603) 


0.14 


9 00+0.02 
z - OJ — 0.01 


-0 1,+ 002 
u - i,3 -0.03 


1.17(814) 


1959+650 


0094380201 


-, 04+O.33 


1 7s +0.25 
J -' 8 -0.25 


1 o 7 +0.32 
i - 8 '-0.31 


1-2918:? 


1.45(57) 


0.31 


1.64 ±0.10 


-0 14+ - 18 
u - i4 -0.19 


1.43(57) 




0094383301 


1 7 cr + 0.19 
J" '°-0.19 


2 14+ ' 12 
z - 14 -0.12 


9 91 +0.09 
z - zl -0.09 


2 S6 +0 16 
z -°°-0.16 


1.11(331) 


0.57 


1 »s+ 004 

1-S8_0 .03 


-0.39 ±0.05 


1.14(513) 




0094383501 


l.86l°;il 


1 oq+0.07 


9 n+0.06 
z - ±1 -0.06 


2.3618:1 


1.05(653) 


0.59 


1.81 ±0.02 


n o 7 +0.03 


1.13(862) 



"'Nh assumed fixed at Galactic (see Table 3) for all fits. 

b X 2 /dof. figures were derived by separately adding the x 2 and dof figures and only then dividing. 
c Four-band not performed because only PN data of low S/N were available 



TABLE 5 

Discrete features found in RGS Spectra 



1ES Name" 


OBSID 


Energy (keV) 


Ax 2 


Width (cV) 


Comments 


1ES0033+595 
1ES1101-232 
1ES1101-232 
1ES1255+244 
1ES1255+244 
1ES1959+650 


0094381301 
0094380601 
0094380601 
0094383201 
0094383201 
0094383301 


1.07 
0.67 
1.36 
1.29 
2.91 
1.92 


-12.02 
-18.56 
-16.62 
-11.98 
-8.91 
-11.77 


45 +20 

40 -24 
7 +4 
'-3 
22 ±9 
30 ± 16 

8±« 
< 18 


MOS 2 + 1 RGS only 
RGS 2 only 
RGS 2 only 

MOS 1, MOS 2 but not in RGS 

Single MOS 1&2 channel, instrumental? 

Two RGS 1 channels, not RGS 2 or MOS 


a This table lists only spectra 


where discrete features 


were found. The spectra 


where the search was 


performed but no features were found 


are listed in §3.3 



TABLE 6 

MULTIWAVEBAND PARAMETERS 



1ES Name 


OBSID 


a RO 


«M 


log Vpea, 


nnQQ i croc; 


UUy4ooloUl 


0.70 


0.66 


17.17 


ni on i qac\ 
U1ZU+D4U 


nnri/i qqoi m 
UUy4oozlUl 


U.00 


u.y / 


"17 "17 


ni A C I 1 QO 

U14o+loo 


nnn /i QO/ini 
UUy4ooo4Ul 


0.59 


1 A A 

1.44 


14.48 


nooQ i noo 
UoZo+UZz 


nnn /i Qooc;ni 


0.49 


1.29 


15. 12 


HQ/I 7 1 Ol 

Uo4 ( — Izl 


nnri/i qqi 1 m 
UUy4oollUl 


n on 


u.yo 


I 7 on 

I I .yy 


ri/i 1 a i nnn 
U414+UUy 


nnri/i qoqi ni 
UUy4ooolUl 


U.oy 


n on 

u.yy 


1 ft QQ 

id. oo 


0o47+250 


0094380901 


0.36 


1.18 


16.28 


0o4 ( +250 


0094382901 


0.36 


1.18 


16.28 


1028+511 


0094381801 


0.42 


0.96 


16.51 




0094382701 


0.47 


1.07 


16.21 


1101-232 


0094380601 


0.45 


0.67 


17.57 


1133+704 


0094170301 


0.26 


1.38 


15.13 


1255+244 


0094383001 


0.53 


0.82 


16.60 




0094383201 


0.53 


0.82 


18.00 


1553+113 


0094380801 


0.32 


1.16 


15.70 


1959+650 


0094380201 


0.47 


0.63 


18.00 




0094383301 


0.48 


0.72 


18.00 




0094373501 


0.48 


0.79 


18.00 



Comments 



Converges poorly; radio flux had to be lowered artificially 



Confusing radio source? Catalog flux 25% of UMRAO value 
a x < 1 - use log v peak =18. 

Catalog optical flux 15x brighter than OM data 
Catalog optical flux 15x brighter than OM data 



Catalog opt flux 25x brighter than OM; increased opt flux 2x 
a x < 1 - use log i/ peak = 18. 

a x < 1 - use log v peak = 18. 

Used log y pea k = 18 even though object not in flare. 
Used log Vp e ak = 18 even though object not in flare. 



